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1 A historical introduction 

Diplon, deuton, deuteron: under different names, the nucleus of deuterium, or 
diplogen, has been the subject of intense studies since its discovery in 1932. 
As the only two-nucleon bound state, its properties have continuously been 
viewed as important in nuclear theory as the hydrogen atom is in atomic theory. 
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Yet, ambiguities remain in the relativistic description of this system and the 
two-nucleon picture is incomplete: meson exchange and nucleon excitation into 
resonances should be considered in the deuteron description. The question of 
rare configurations where the two nucleons overlap and loose their identity is 
still under debate. We are still looking for the elusive effects of quarks in the 
nuclear structure. 

In this year of the millenium, the present article will first attempt to recall 
the early discoveries, measurements and theories. It will then boldly jump over 
decades of continuous efforts, building upon these, to present not an exhaustive 
review but an up-to-date status of our understanding of the deuteron, with a 
special emphasis on its electromagnetic form factors. To do justice to some 
seventy years of activity in this field is an immense task which is more easily 
approached by quoting here several reviews along this path |[ ^, |4[ ^, ||, 0, 
|[ ^ |l^, |l^, |l^, The important subject of electro- and photo-disintegration 
of the deuteron will be only partly covered, referring the reader to ||lj] . 



1.1 Discovery of the deuteron 

The existence of the first isotope of hydrogen was suggested in 1931 by Birge and 
Menzel lls] in order to remove discrepancies between two different measurements 
of the atomic mass of hydrogen. A first estimate of an abundance ratio ^H/^H = 
4500 was inferred from this hypothesis, close indeed to the actual value of 6700. 
The stable isotope was discovered by Urey and collaborators |l6j a few months 
later, investigating distilled samples of natural hydrogen for the optical atomic 
spectrum of in a discharge tube. Isotopic separation to study the properties 
of deuterium quickly became an intense activity. Its mass was measured by 
Bainbridge [0. While Chadwick was discovering the neutron, several tens of 
papers were written, in one years time, devoted to the study of deuterium. 
An illuminating summary of this early research was made by Bleakney and 
Gould |18[. In 1933, "deutons" were used as accelerated projectiles first at 
Berkeley |]19t, then at Caltech and at Cavendish. Chadwick and Goldhaber |^ 
measured the first photodisintegrations pn in 1934. 



1.2 Early theories 

In 1932, there was no satisfactory theory of the nucleus. The nucleus was 
thought to be composed of protons and electrons since these were the only 
known charged particles and nuclei were seen to emit electrons (/3 decay). The 
electrons were needed to cancel the positive charge of some of the protons in 
order to account for nuclei with identical charges, but with different masses, and 
to allow for the possibility of binding of the nucleus by means of electric forces. 
This was clearly unsatisfactory because the Coulomb force could not account 
for the binding energies of nuclei and the attempt to construct the nuclei from 
the incorrect number of spin- 1/2 particles could not produce the correct nuclear 
spins. 
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The discovery of the neutron, shortly after that of the deuteron, did not im- 
mediately eliminate the confusion since the previous model persisted by simply 
describing the neutron as a bound system of a proton and an electron. Based on 
this faulty assumption, Heisenberg produced the first model of proton-neutron 
force |2^. Since it was not possible to actually construct a description of the 
neutron with the ep model, Heisenberg simply assumed that the pn force could 
be described by a phenomenological potential and that the neutron was a spin- 
1/2 object like the proton. Based on an analogy with the binding of the 
ion by electron sharing, Heisenberg proposed that the force must involve the 
exchange of both spin and charge in the form of Forces con- 

taining the remaining forms of spin and isospin operators were soon introduced 
by Wigner Majorana |Q and Bartlett In all cases the spatial form 

of the potentials was to be determined phenomenologically to reproduce the 
deuteron properties and the available nucleon-nucleon (NN) scattering data. 
In 1935Bethe and Peierls wrote the Hamiltonian of the "diplon" with an 
explicit introduction of a short range interaction. This approach became the 
mainstay of nuclear physics which has produced considerable success in describ- 
ing nuclear systems and reactions. The ep model of the neutron was not com- 
pletely abandoned until after the Fermi theory of (3 decay became widely 
accepted. 

The progress in discoveries and understanding was then so great that, in spite 
of an otherwise bleak social or political situation in many countries involved, this 
period is recalled as "The Happy Thirties" from a physicist's point of view pTf . 

One of the other great theoretical preoccupations of the late 1920's and the 
1930's was the development of quantum field theory starting with the first works 
of Dirac on quantum electrodynamics (QED) ||2^, the Dirac equation for the 
electron and the Dirac hole theory |3^ with field theory reaching its final 
modern form with Heisenberg |3l| . QED at this time was very successful at 
tree-level but the calculation of finite results from loops was not really tractable 
until the introduction of systematic renormalization schemes in the late 1940's. 
The first attempt to apply quantum field theory to the strong nuclear force was 
Yukawa's suggestion |^ that the force was mediated by a new strongly cou- 
pling massive particle which became known as the pion. This started another 
strong thread in the theoretical approach of the nucleus by using meson-nucleon 
theory to obtain nuclear forces consistent with the phenomenological potential 
approach. The primary attraction of this approach is that a more microscopic 
description of the degrees of freedom of the problem is provided and that addi- 
tional constraints are imposed on the theory by the necessity of simutaneously 
describing nucleon-nucleon and meson-nucleon scattering. Ultimately, as it be- 
came clear that the mesons and nucleons were themselves composite particles, 
meson-nucleon theories were replaced as fundamental field theories of the strong 
interactions by quantum chromodynamics (QCD). However, the meson-nucleon 
approach is still a strong element in nuclear physics as a basis for phenomenol- 
ogy and is making a potentially more rigorous comeback in the form of the 
effective field theories associated with chiral perturbation theory. This situa- 
tion is unlikely to change until it becomes possible to at least describe the NN 
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force and the deuteron directly from QCD. 



1.3 Spin 

Breit and Rabi |Q first suggested the use of magnetic deflection of an atomic 
beam in an inhomogeneous field to measure nuclear spins. The coupling of 
electronic (Je) and nuclear (J) spins is not totally negligible compared to the 
coupling of the electronic spin to the external magnetic field, provided the latter 
is weak enough. One then observes (2 Jg + 1) x (2 J+ 1) lines with a predicted in- 
tensity pattern. The atomic and molecular beam studies were to be implemented 
with great success (see Sec. ||), but the first determination of the deuteron spin 
used other methods. 

Farkas and collaborators |Q demonstrated the ortho-para conversion in the 
diplogen (as they called the deuterium molecule) and determined the spin and 
statistics of the nucleus from the equilibrium ratio between these two states at 
different temperatures. They concluded that the diplogen nucleus must obey 
Bose- Einstein statistics, that the most probable value of its spin was 1, and that 
its magnetic moment was about one fifth of that of the proton. 

Using photographic photometry, the alternating intensities in the molecular 
spectrum of deuterium were investigated by Murphy and Johnston psf , who 
concluded that indeed J = 1 for the "dcuton" . 



1.4 Connection with OPE 

The deuteron thus quickly appeared as a loosely bound pair of nucleons with 
spins ali gned ( spin triplet state). The existence of a small quadrupole moment 
(see Sec. 3.1.3 ) implies that these two nucleons are not in a pure S state of rela- 
tive orbital angular momentum, and that the force between them is not central. 
Taking into account total spin and parity, an additional D wave component is 
allowed. Such a D wave can be generated by the tensor part of the one-pion 
exchange (OPE) potential [|, M. 



2 The nonrelativistic two-nucleon bound state 

2.1 The potential model of the deuteron 

The potential model of the deuteron is described by the Hamiltonian 

H = fi+f2 + V (1) 

where Ti is the kinetic energy operator for particle i and V is the two-body 
potential. Successful NN potentials must have several basic characteristics in 
order to satisfactorily describe the deuteron static properties and the NN scat- 
tering data. The long distance part of the potentials is described by one-pion 
exchange while the intermediate and short range parts may be either param- 
eterized in terms of simple functional forms, or obtained from models involv- 
ing meson exchanges. The very strong anticorrelation of nucleons requires that 
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these potentials be repulsive at short distances. The potentials must have terms 
involving scalar, spin-spin, tensor and spin-orbit forces. The tensor force is of 
particular importance in producing the single spin-1, iso-singlet deuteron bound 
state. The long range tensor force is provided automatically by the exchange 
of the pseudoscalar pion. Modern phenomenological potentials also include ad- 
ditional nonlocalities by means of terms quadratic in the relative momentum 
and/or quadratic spin-orbit terms. Improved fits to scattering data also require 
that isospin symmetry breaking be imposed via the inclusion of electromag- 
netic interactions between nucleons and by additional explicit isospin symmetry 
breaking terms in the potential. By fitting the potentials directly to the scat- 
tering data, several phenomenological potentials have been contructed that fit 
the scattering database with very close to I. 

An discussion of the most commonly used NN potentials may be found in 
the review These include the so-called Reid-SC |3^, Paris Bonn |]39| , 
CD-Bonn Nijmegen Reid93 Q and Argonne vis @ potentials. 

Given a potential, the resolution of the Schrodinger equation in the T = 
0, J = I np channel leads to the bound state wave function discussed hereafter. 



2.2 The deuteron wave function 

The tensor force requires that the deuteron wave function be a mixture of '^Si 
and '^Di components, so the deuteron wave function is of the form 

(x) = !^yM^ ie,<^) + ^y,'i, {9,^), (2) 

where 

yjlsieA)= {J,M\L,mL;S,ms)YLM{e,cly)\S,ms) (3) 

are the spin-spherical harmonics. The reduced radial wave functions u(r) and 
w{r) correspond to the S and D waves respectively. The S and D state proba- 
bility densities are defined as 

pg(r)=u^(r) and poir) = (r) . (4) 

The corresponding S and D state probabilities are then 

poo POO 

Ps = ps{r)dr and Pd = PD{r)dT (5) 
JO Jo 

and the normalization of the wave function requires that 

Ps + Pd = 1. (6) 

The reduced radial wave functions for the Argonne vis potential are shown in 
Fig. |. The wave functions for other modern potentials are very similar. 
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Figure 2: The deuteron S wave function in configuration space and in momen- 
tum space: u{r)/r and pu{p) (calculated from the Argonne wis potential). 



From the wave function, a characteristic size of the deuteron r„i is defined 
as the rms-half distance between the two nucleons : 



/•oo 







-u^(r) + w^(r)] r^dr (7) 



A conspicuous feature of the nuclcon-nucleon interaction is the short range 
repulsion, which leads the radial S wave function u{r)/r to be significantly re- 
duced at distances smaller than approximately 1 fm (see Fig. ||). This introduces 
a distance scale in the wave function in addition to the overall deuteron size. 
This small distance behaviour is the subject of most of the experimental and 
theoretical studies which will be presented in Sees. || and |[ As a result of this 
dip at small r, the Fourier transform u{p) contains a node at approximately 2 
fm~^, as seen also in Fig. ^j. The u and w wave functions are given in momentum 
space by : 

/•oo />oo 

"(p) — / u{r)jo{pr)rdr and w{p) = — / w{r)j2{pr)rdr . (8) 
Jo Jo 

To conclude this presentation of the size and shape of the deuteron, the densities 
[^"(x)!^ and |V'^(x)p are illustrated in Fig. ||. 
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Figure 3: Deuteron densities in M = (left) and M ~ 1 (right) magnetic 
substates. The red spots correspond to the maximal nucleonic densities, while 
the dark volumes correspond to lower densities (outer surface is for 10% of 
maximal density). See for equivalent representations. 



In this simple potential model of the deuteron, the magnetic moment of the 
deuteron is determined entirely by the D state probability Pd ■ 



fJ-d = fJ-s 



2 IM. 



(9) 



where fig = fJ-n + fJ'p is the isoscalar nucleon magnetic moment. The deuteron 
electric quadrupole moment is also determined from the wave functions: 



Qd 



1 



w{r) 



u{r) 



1 

71 



w{r) 



r dr . 



(10) 



In both cases, these quantities are modified by extensions to the basic potential 
model (see Sec. H). In particular the direct relationship between the magnetic 
moment and the D state probability will be broken by such extensions and, 
therefore, this probability is not an observable Q. 

Since the nuclear force is of finite range, it is easy to determine the asymp- 
totic form of the wave functions 



u{r) ^ Ase ''^ and w{r) ^ A^e 



■'yr 



1 



3 



dry 



(11) 



where 7 ~ y^em, with m being the reduced np mass and e the deuteron binding 
energy (see Ref. 45 for a relativistic definition of 7). As and Ad are the 
asymptotic normalization factors, determined by matching the asymptotic form 
( pj] ) to the calculated wave functions in the interior region where the potential 
is nonvanishing. As and the ratio 



Vd 



Ad 
As 



(12) 



are directly related to observables as discussed in the next section. 
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Table 1: Experimental determinations of the deuteron static properties 



Quantity 


Most recent 
determination 


Value 


Mass Md 


|4 


7, 


48 


4875.612762 (75) MeV 


Binding energy e 




49 




2.22456642 (48) MeV 


Magnetic dipole moment fid 




48 




0.8574382284 (94) 


Electric quadrupole moment Qd 


[46, 


50 


. 51| 


0.2859 (3) fm^ 


Asymptotic ratio rjd = Ao /As 




52 




0.0256 (4) 


Charge radius rch 
Matter radius 




53 






2.430 (40) fm 
4.975 (3) fm 


|5 


4. 


55 




Electric polarizability as 


|5 


5, 


57 




0.645 (54) fm^ 



3 Static and low energy properties 

A review of the measured static properties of the deuteron, together with the 
low energy neutron proton (np) scattering parameters, was given by Ericson 
and Rosa-Clot |6[ who studied in detail their connection with NN potential 
models. An updated status of the experimental information on the deuteron is 
given in Table |4| and discussed hereafter. 

3.1 Deuteron static properties (experiment) 



3.1.1 Mass and binding energy 

In the past ten years, significant progress in the precision of the measurement 
of some atomic masses has been made by comparing cyclotron frequencies of 
different pairs of ions in a Penning trap. In this way, the deuterium atomic mass 



is known with a relative precision of 40^ 47 1. 

The deuteron binding energy is best determined by measuring the energy of 
the gamma-rays coming from radiative np capture with thermal neutrons. This 
energy is now measured to a relative accuracy of 2 x 40~^, using a crystal diffrac- 
tion spectrometer At this precision, even for such a low energy process, 
some of the earlier work may have to be corrected for relativistic kinematics and 
the Doppler effect 

The binding energy and the mass measurements can be combined for the 
most precise determination of the neutron mass |^9|. The errors on the val- 



ues reported in Table |4 
934.494043(37) MeV { 



include the uncertainty in the atomic mass unit u = 
. They are well beyond the accuracy of nuclear models. 



3.1.2 Magnetic dipole moment 

The first measurement of the deuteron magnetic moment was performed by Rabi 
in 1934 fs^], based on a principle already alluded to. From the deflection 
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of an atomic beam in an inhomogeneous magnetic field to the use of molec- 
ular beam resonance and other methods, these techniques were continuously 
improved [ |60| . Precise measurements of nuclear magnetic resonance frequencies 
of the deuteron and proton in the HD molecule give the ratio of deuteron to 
proton magnetic moments. However, the adopted value in Table |^ results from 
a simultaneous determination of the electronic and nuclear Zeeman energy lev- 
els splittings in the deuterium atom, yielding the ratio of deuteron to electron 
magnetic moments [|8|. 

3.1.3 Electric quadrupole moment 

The deuteron was found to possess an electric quadrupole moment in 1939 | ]6l[ . 
This discovery had far reaching consequences: it meant that nuclear forces were 
not central and were more complex that previously thought. It was to become 
the best qualitative and quantitative evidence for the role of pions in nuclear 
physics 

In contrast to the case of the magnetic moment which is determined through 
its coupling to an external applied magnetic field, the quadrupole moment does 
not couple to an external electric field. One measures instead, in HD or D2 
molecules, the interaction of the deuteron quadrupole moment with the elec- 
tric field gradient created along the molecular axis by the neighbouring atom. 
The experiment provides an electric quadrupole interaction constant which 
must be divided by the theoretically calculated field gradient |0 to obtain the 
quadrupole moment. 



3.1.4 Asymptotic ratio D/S 

The ratio (^|) is deduced from measurements of tensor analyzing powers in 
sub-Coulomb reactions on heavy nuclei by comparison with calculations 

in the distorted wave Born approximation (DWB A) |52f| . The value of rjd is then 
directly proportional to the analyzing powers. Other determinations based on 
dp elastic scattering rely on pole-extrapolation and are somewhat less precise. 



3.1.5 Radius and size 

The deuteron size may be characterized by a charge radius r^-h and by a matter 
radius rm- The latter is defined from the deuteron wave function ( [^). 

Elastic electron scattering has been used since the early fifties ^21 to mea- 
sure the shape of nuclei. This topic will be discussed at length in Sec~§. At low 
momentum transfer, the cross section data yield the charge rms-radius of the 



target nucleus through the relation r^^ = —6 rfCc/c^Q^lna^n (^^® ^JJ ^^"^ 



' ch — uvj^GV"^ |q2^0 

the definition of Gc)- It was demonstrated recently that a precision extraction 
of the deuteron rms charge radius from electron scattering data requires taking 
into account the Coulomb distortion of the incoming and outgoing electrons |^ . 
A new analysis of the world data was then performed, yielding the value in Ta- 
ble ^ The quoted uncertainty combines quadratically the fit statistical error 
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and the dominant systematic error, the latter coming mostly from experimental 
normalization uncertainties. The usual radiative corrections to electron scatter- 
ing do not include the contribution of hadronic vacuum polarization, but that 
effect should be smaller than the present uncertainties when extracting charge 
radii Q. 

The rms charge radius Vch is related to the matter (or rather nucleonic) 
rms-radius by 

''c/. - + Ar2„ + r2 + r2 + ^ ^ (13) 

where rp = 0.862(12) fm is the proton charge rms-radius |^ and = —0.113(5) 
fm^ is the neutron charge ms-radius p6t . Ar^ is a contribution from non- 
nucleonic degrees of freedom, close to but with an uncertainty estimated to 
±0.01 fm^. The quantity given by = r^„ + Ar^ is usually defined as the 
deuteron radius. The last term in ( p^ ) is of relativistic origin Note that 
the above quoted value of rp, as extracted from ep elastic scattering, is in slight 
disagreement (2cr difference) with recent Lamb shift measurements Finally, 
the theoretical uncertainty in the deuteron radius associated with the correction 
due to the nucleon finite size has been estimated to about 0.002 fm ]69[ |. 

The nuclear-dependent correction to the Lamb shift in hydrogen and deu- 
terium atoms is directly proportional to the nuclear mean-square radius. From 
the isotope shifts in the pure optical frequency of 15 — 25* two-photon transitions 
in atomic hydrogen and deuterium, the difference r^^ — of mean-square charge 
radii for the deuteron and proton is accurately determined p^ . Small correc- 
tions due to the deuteron polarizability seem to be under control [55[ |. Then 
from (|l^ ) and the value of r^^, is extracted with a better precision than Tch 
from ed scattering. Our quoted uncertainty is larger than in Ref. because 
of the use of a larger uncertainty in and the addition of the uncertainty due 
to Ar^. ^ ^ ^ ^ _ 

Taking into account additional small corrections summarized in Ref. [p3[ , the 
two results given in Table are quite compatible, in the sense that they satisfy 
Eq.(^|). Furthermore, the value of r„i follows the expectations from modern 
NN potentials. 

3.1.6 Electric polarizability 

The electric polarizability ue characterizes how the deuteron charge distribution 
can be stretched and acquire an electric dipole moment under the influence of an 
external electric field. It was determined through elastic scattering of deuterons 
from "^^^Pb well below the Coulomb barrier and extracted from low energy 
photoabsorption p7|. The two results are slightly incompatible {2a difference). 
The value in Table |l| is our average. 
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3.2 Low energy np scattering parameters 

The deuteron may also be viewed as a pole in the S'-matrix describing the 
np scattering in the coupled ^Si and ^Di channels. This S'-matrix can be 
experimentally determined from a phase-shift analysis of the scattering data. 
An extrapolation to negative energies down to the measured deuteron binding 
energy, either by an effective range expansion |^ or a P-matrix approach jTOf , 
allows to extract rjd and As ■ The asymptotic ratio rjd — 0.0254(2) is given by the 
extrapolated mixing parameter ei while the asymptotic S state normalization 
As = 0.8847(8) fm~^/^ is essentially related to the effective range and thus to 
the triplet scattering length at (numerical values from |7l|). 



3.3 Static properties and the A^A^ potential 

All deuteron static properties discussed above are well reproduced by NN po- 
tential calculations such as Argonne wis, Nijmegen II, Reid93 or CD-Bonn, with 
the notorious exception of the quadrupole moment, which is always a few per- 



cent too low (see also Sec. 5.1). Meson exchange contributions, to be discussed 
later, must be taken into account for a better agreement with data. The binding 
energy e is taken as a constraint in the determination of all potentials. 

Compilations of deuteron static properties caculated with recent NN inter- 
action models appear in Note that most potentials result in a D wave 
probability Pd between 5.6 and 5.8%, except for the CD-Bonn potential where 
Pd = 4.83%. 

Various correlations were established between the calculated static proper- 
ties, independently of the NN potential used. For example, linear relationships 
between As and or A|(l-|-77^) and r^j |5|], and Qd/A'g and rj ^ were es- 

tablished, the latter depending on the value of the ttNN coupling constant used 
in the potential calculation. For more recent potentials, linear relationships be- 
tween /id and rj on one hand, Qd and r/ on the other hand, are illustrated in [fz^ . 
Finally the electric polarizability a^; is directly proportional to |55|]. 



4 Elastic electron-deuteron scattering 

4.1 Deuteron electromagnetic form factors 

Invoking Lorentz invariance, current conservation, parity and time-reversal in- 
variance, the general form of the electromagnetic current matrix element for 
elastic electron scattering from the spin-1 deuteron can be shown to have the 
general form j7|] : 



G^,,_^(P',P) = -|Gi(Q^)(ei,(P')-a.(^))(/^' + n^ 
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- Gs{Q')^jCx'^{p')-qmAP)-q)ip' + prj (14) 

where Md is the deuteron mass, P and P' are the initial and final deuteron 
four- moment a, q = P' — P is the virtual photon four-momentum, (P) and 
^x*{P') ^I'G the polarization four- vectors for the inital and final deuteron states. 
The Gi{Q^) are form factors depending only upon the virtual photon four- 
momentum; assuming hermiticity, they are real. Since the virtual photon four- 
momentum is always spacelike for electron scattering, we use the convention 

The current may be expressed in terms of charge monopole, magnetic dipole 
and charge quadrupole form factors. These are related to the Gi{Q^ys by: 



GciQ') = GiiQ^) + -vGQiQ') 

Gm(O') = G2(g2) 
Gq(Q') = G,{Q^) - G2{Q') 



ii + v)G3iQ' 



(15) 



with 



These form factors are normalized such that 



Gc(0) - 1, 

Ga/(0) = ^fid: 



Gq(0) 



M^Qd 



(16) 



(17) 



The experimental values of Gj\/(0) and Gq(0) are respectively 1.714 and 25.83 
(see Table 0). 

4.2 Observables 

In the Born approximation of a onc-photon exchange mechanism and neglect- 
ing the electron mass, the cross section for elastic scattering of longitudinally 
polarized electrons from a polarized deuteron target can be calculated from the 
current to give in the laboratory frame |74[| : 



da 

dn 



^ sin^ ^ 

Md 2 



[vlRl + vtRt + vttRtt + vtlRtl 

+ 2hvT'RT' + 2hvTL'RTL'] , (18) 



where 



a cos -, 
2E sin^ 



(19) 
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is the A/[ott cross section, E the electron beam energy, 9 the electron scattering 
angle and h = ±^ the electron helicity. The Ri are response functions and the 
Vl kinematical factors are: 



For elastic scattering 



Vl = 






Vt = - 


+tan - 


vtt = 


2 q2 




vtl = - 




q2 


Vt' = 


— + tan - 
q^ 2 


VtL' = - 


1 Q\ 






1 




q2 1 





+ tan^ 



tan ■ 



Each of the response functions can be written as 

9 



i=l 



(20) 



(21) 



(22) 



where I = {L, T, TT, LT, T', LT'} and the n are members of the set of the 
unique deuteron density matrix elements expressed in terms of elements of a 
spherical tensor pkq- This set is represented by 

n = |poo,y^Pio,-^P20,\/3-Rep22,\/3-fmp22, 



-Repi-i, \ -Impii, \ -Rep2i, \ -Imp; 



'J21 



(23) 



The nonvanishing reponse functions for elastic scattering may be written in 
function of the deuteron form factors : 



Rl{Q^Poo) 



V2 

Rt{Q\poo) 



vGciQ')GQ{Q') + -v'GliQ') 



-rjil+r,)GUQ') 
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--rj{l + rj)GUQ') 



22) 



Rt'{Q^ 



Pw) 



RTL'iQ,\-Repii) = 



-^775 (1 + rj)i Gm{Q^) \Gc{Q^) + ^Gq{Q%A) 



da 



It is conventional to write the cross section as 



1 + 2B - 2 9 



X 



A(Q^)+S(g^)tan^ 



X [1 + P20 ■ i2o(Q', 0) + 2i?ep2i • t2i(Q', e) + 2Rep22 ■ t22(Q', 0) 
+hpw ■ hoiQ^ 0) + 2hRepn ■ tn{Q^ 0)] (25) 

where the unpolarized elastic structure functions A{Q'^) and B{Q^) are defined 
as 



A(Q2) = vlRl{Q\U) 



1 



2(1 + 77) 



i?T(Q',C/) 



= Gl{Q^) + tr,Gl,{Q') + %'Gl{Q') 



(26) 



and 



B{Q') = Rt{Q,U) 



-r,(l + rj)Gl,{Q^). 



(27) 



Note that the dependence of the cross section on the target polarization is 
conventionally given by analyzing powers denoted T^g; we implicitly use here 
the equivalence between analyzing powers and recoil deuteron polarizations: 
tkq = Tkq . Defining 

Ro = A{Q^) + B{Q^)tan^^, 



(28) 



the tensor polarization observables are 

1 



t2o{Q,0) 



\/2i?o 

{j^vGc{Q')Gq{Q') + Iri^GliQ^) 



1 + 2(1 + 7/) tan^ 



(29) 
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ZUo 



2V3R0 



(30) 



t2i(g' 



V^Ro'^i ^ + ^(1 + V) tan2 - GmIQ") Gq(Q2) (31) 



tio(Q',e) 



-^/|/(l + .)tany^+tan^^GL(Q^) (32) 



tii(Q',e) = 







VSRo 



V^H^ tan -Gm(Q') Gc(Q') + :^Gq(Q2) (33) 



There are only two unpolarized elastic structure functions, but three form 
factors. Since B{Q^) depends only on G\j[Q'^), this form factor can be deter- 
mined by a Rosenbluth separation of v4(Q^) and B{Q^), or by a cross section 
measurement at 9 = 180°. A{Q^) depends on all three form factors so that only 
a quadratic combination of Gc{Q^) and Gq{Q^) (the longitudinal part of A) 
can be determined from the unpolarized cross section. A complete separation 
of the form factors therefore requires the measurement of at least one tensor 
polarization observable. The possible candidates are i20j ^21 and tu, since i22 
and tio depend only upon G\,j{Q'^) and the unpolarized structure functions. ^21 
and ill, being both proportional to Gm(Q^), are of smaller magnitude than 
^20 and provide in practice a smaller "lever arm" to determine either Gc{Q^) 
or Gq{Q'^) In addition, the measurement of tu requires intense polarized 
electron beams, which became available only recently, but offers the simplifica- 
tion of a vector polarization measurement. In all cases so far, this leaves ^20 as 
the observable of choice to extract Gc{Q^) and Gq{Q^). 

Note that all of the polarization observables depend upon the scattering an- 
gle 9 through kinematical factors and Ro{Q^,9). Consequently, the polarization 
observables measured in different experiments under different kinematical con- 
ditions can only be compared if some convention is assumed. Since the first 
t2o measurement |Q was performed close to 9 = 70°, it has been customary 
to quote the observables at that angle. For experiments not performed at 70°, 
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the observable is extrapolated to this angle, using the known A{Q'^) and B{Q^). 
Another convention is to use the alternate quantity: 



<20(0") 



2\ — 



1 i?i(g2,-l=T2o) 



V2 Rl{Q^U) 



V2 



Gim + WGliQ^) 



(34) 



The choice of f20i which is not strictly speaking an observable, has several 
advantages. t2o is independent of 9, and thus depends only on Q^. It is a purely 
longitudinal quantity, and as such is independent of the magnetic form factor. 
It has simple properties related to those of the charge and quadrupole form 
factors (see Sec. iA and Refs. j?^, |7^). In particular, the position of nodes in 
these form factors may be determined directly from a plot of i20- Numerically, 
i2o{Q^) can be determined from A(Q^), B{Q'^) and t2a{Q^,9) through 



i2o(Q') 



^2o(g^e) + 5^ 



1 



with S ■ 



vt 



(35) 



As illustrated in Sec. [4.3| , in all measurements to date, the ratio S is small, so 
that t2o{Q^), t2o{Q^, 9) and t2o{Q^, 70°) are not very different from each other. 
For all these reasons we will use this quantity along with A{Q'^) and B{Q^) for 
comparison of theoretical predictions to data. 

In closing this introduction to ed elastic scattering observables, we refer to 
App. A for a short discussion of a possible two-photon exchange contribution, 
especially in view of the large range of available data. 



4.3 Review of elastic ed data 



The first experiment to measure elastic scattering of electrons on the deuterium 
was performed at the Stanford Mark III accelerator jVOj. Since then, many 
cross section data points have been measured at various accelerators over the 
world, with ever increasing precision and at larger and larger momentum trans- 

fers |§, m, n, n, ||, |], g, n g|, n, g, n HI], n, 

Quite spectacular is the recent achievement of Jefferson Lab to measure A(Q^) 
up to ~ 6 (GeV/c)^, making use of a record luminosity of about 5 x 10'^^ 
cm~^s~"'^ to reach cross sections as low as 10"'*^ cm^/sr New measurements 
of B{Q^) for ~ 0.7 to 1.3 (GeV/c)^ from the same experiment will soon be 
available ||lO0|] . The kinematics of all these experiments are illustrated in Fig. ^. 
Forward angle scattering yields the elastic structure function A, while back- 
ward angle scattering allows the determination of the elastic structure function 
B. The dashed lines in Fig. ^ indicate what fraction of the cross section cor- 
responds to the contributi on o f B{Q ^) tan'^(6'/2). Other experiments measured 
cross section ratios ed/ep [101, 102|. 
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Figure 4: Kinematical settings [0 vs Q) for various cross section measure- 
ments in ed elastic scattering. The small angle data from SLAC |90) (open 
circles) and JLab/HallA (filled circles) extend respectively to 10.1 and 12.4 
fm"-"^, beyond the scale of the horizontal axis. The dashed lines correspond to 
Btan2(6'/2)/i?o =0.01, 0.1, 0.5, 0.9, 0.99, respectively from bottom to top. 
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As already mentionned, the separate determination of the deuteron charge 
monopole and quadrupole form factors necessitates the measurement of a polar- 
ization observable. The observable of choice is t2o, which is a measure of the rel- 
ative probabilities of finding the deuteron in magnetic substates M = — 1, 0, -|-1 
after the ed scattering. The Bates Linear Accelerator Center was the first to 
measure ^20 | |76| and to provide an experimental evidence for the existence of 
a node of the charge form factor |7^. These double scattering experiments 
were recently brought to the (up to now) highest possible momentum transfers 
at Jefferson Lab [103|. For this last experiment as well as for the previously 
mentioned A measurement |Q, electron beams of 100 to 120 /xA were used in 
conjunction with liquid deuterium targets up to 15 cm long, capable of dissipat- 
ing 600 W of power deposited by the beam. The combination of a record inte- 
grated luminosity, in excess of 10^ pbarn"^, and of a large acceptance magnetic 
channel focusing recoil deuterons onto the high efficiency polarimeter POLDER 
(see Fig. ^ and App. B) allowed the measurements of t2o to be extended up 
to ~ 1.7 (GeV/c)^. The alternative measurement of the tensor analyzing 
power T20, using a polarized target intercepting a stored electron beam, was ini- 
tiated at Novosibirsk |]lO|, |l0|, |l06), and improved at NIKHEF flO^, |l08). New 
preliminary results from Novosibirsk |10!;] are also included in this review. On 
the other hand, the use of a solid cryogenic target (ND3 ) in an external electron 
beam at Bonn resulted in a too low luminosity | [110| . The kinematical settings 
of all these experiments are illustrated in Fig. Id. In all cases, the magnetic 
contribution to i20 is small. A further comparison between these polarization 
measurements is contained in App. B. 

The other tensor polarization observables t2i and t22 (or T22) were also 
measured ^ ^ ^ ^, ^ 

The A data at low and high 
^, in particular in Figs. ^ 
The ^21 data appears in Fig. 27. In closing this section, let us mention 
a few inconsistencies in this data set in the light of recent measurements. [| The 
^(Q^) Cambridge |^ and Bonn Q data are very probably too low, since both 
recent measurements at Jefferson Lab Q agree with the "higher" trend al- 
ready given by the SLAC data |90|. Still, as apparent in Fig. these two 
JLab measurements differ from each other (10-15%) in the region = 1 to 2 
(Ge V/c) ^. For a comparative discussion of these two independent experiments, 
see |lll| . The t2o (or T20) data are necessarily less precise than cross section 
measurements and naturally exhibit some scatter. Although all data points 
are compatible with parameterizations such as discussed below, there are some 
trends between different data sets. No parameterization or model can accomo- 
date both the Bates ||7^ and NIKHEF |108| data sets: one or the other is too 
low, or both arc. The same Bates data are also systematically lower than both 
the JLab |103| and the preliminary Novosibirsk data |109|, and the precise low 




Figure [7| shows a good part of the existing data 
will be better illustrated in the figures of Sec 

and 



^ At the time of print of this paper, a better determination of the beam energies at 
JLab/Hall C results in some corrections, within quoted errors. The A{Q^) values Q should 
decrease by 1 to 3% with increasing Q^, while the first t2o point [103| should move down by 
about one third of its error. These numbers are subject to confirmation. 
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Figure 5: Experimental set-up for the recent i20 double scattering experi- 
ment |103| in the Hall C (50 m diameter) of Jefferson Lab: e + ^(LDj) 
e' (HMS) -I- d (Deuteron channel) followed hy d + p {pp)n in the POLDER 
polarimeter. 



NIKHEF point |107| is lower than many theoretical expectations. Though 
these scatters are compatible with the quoted experimental errors, they demon- 
strate (together with the theoretical models to be discussed) a need for a more 
accurate measurement in the region Q = 3 to 4.5 fm~^. Finally, the Bates ^22 
data point at Q ~ 3.8 fm^^ |^ is obviously wrong, but the authors did not 
find a correlation between t2o and ^22 in their analysis . 



4.4 Empirical features of form factors 

The three deuteron electromagnetic form factors may be calculated at a fixed 
value of from measurements of A, B and i20- The magnetic form factor Gm 
is readily available from the B measurements, while the two charge form factors 
Gc and Go are determined from vlRl{Q'^, U) = Al ^ A - B/2{1 + 77) <^ 
and i20 ( |34^5| ). The resolution of these equations is most simply described in 
Ref. , together with a discussion of possible ambiguities in the choice of dif- 
ferent solutions in the Q-regions where ^20 reaches its extrema. This procedure 
allows a direct comparison of theoretical models with the form factors, instead 
of observables, but it is limited to the domain where the three observables are 
measured, which is g = - 7 fm-i(Q2 < i.g (GeV/c)^). 

The most striking result is an experimental determination of the node of 
Gc, at Q = 4.21±0.08 fm"i @, which corresponds to ^20= -1/a/2 . This 
behaviour of Gc, though expected from most models, could not have been seen 
with cross section measurements only. The exact location of this node is sensitive 
to the strength of the NN repulsive core (in the impulse approximation, it is 
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Q (fm"') 



Figure 6: Kinematical settings for measurements of elastic ed tensor polarization 
observables. Filled symbols for ^20 measurements of Ba tes/ Argonne f76| (cir- 
cles), Bates/AHEAD (squares) and JLab/POLDER [|03| (diamonds), each 
at a fixed recoil deuteron angle. Open symbols for T20 measurements of Novosi- 
birsk ||l0|, ^ (triangles up), (lO| (triangles down), NIKHEF [|07| (circle), 
[ 108 (squares) and Bonn |11C|| (star), each at a given electron beam energy. The 
tilted bars indicate the detector acceptances. The dashed lines are a measure 
of the small magnetic contribution to t2o'. they correspond to ^20 ~ ^20 = 0.05, 
0.1, 0.15, respectively from bottom to top. 
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Figure 7: Elastic ed scattering observables. The A data are from Cambridge p7| 
(triangles), SLAC Q (open circles), Bonn Q (stars), Saclay |Q (open dia- 
monds), JLab/HallA [|| (full circles) and JLab/HallC @ (full diamonds). 
The B data are from Saclay (diamonds), SLAC [|6) (circles), Bonn |||| 
(stars). For t2Q data legend, see Fig. ||. 
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connected to the node of u{p) in Fig. ^) and to the size of relativistic corrections 
and of the isoscalar meson exchange contributions. A secondary maximum is 
also determined from the data ]103[ |. Thus, hke any other nucleus, the deuteron 
appears to have a charge form factor with an oscillatory diffractive pattern. 
Unlike other nuclei, the sign of this form factor is determined as well. 

The quadrupole form factor Gq exhibits a monotonous exponential fall-off, 
and its first node, corresponding to the yet unobserved second node of ^20: is 
expected beyond Q — 7 fm~^. 

Finally the magnetic dipole form factor Gm has a node at Q — 7. 2zfc0.3 
fm~^, determined by B data and in a lesser extent by t2i data |103|. The 
position of this node is very sensitive to small non-nucleonic components in the 
deuteron wave function, to nucleon-nucleon components of relativistic origin, as 
well as to the description of the pirj contribution to be discussed in the next 
section. 

The three form factors were parameterized in three different ways: rational 
fractions (I), sum of Lorentzian functions in an helicity basis (II) and sum of 
Gaussian functions (III) |78), fitting directly the measured observables, i.e. dif- 
ferential cross sections and polarizations. Figure ^ gives a representation of the 
form factors with an updated version of parameterization I, taking into account 
the preliminary T20 data of Ref . ]109(| . The data base and the parameterizations 



are available in |112|. Another representation of the experimental knowledge of 
the deuteron form factors is given in Fig. ^, where the contributions from each 
of them to the elastic structure function A are given, calculated using parame- 
terization I. 



5 Theoretical issues 

The different classes of theoretical models of the deuteron electromagnetic form 
factors are presented here, illustrated with the most recent calculations on the 
subject. A summary of earlier theoretical work on the subject may be found for 
instance in . In all figures of this section, the data legend is the same as 



in Sec. 4.3 



5.1 Deuteron elastic form factors in the simple potential 
model 

Once the wave functions for a given potential model are obtained as discussed 
in Sec. |^, the electromagnetic form factors may readily be calculated in the 
nonr elativistic impulse approximation (NRIA) using a well established formal- 
ism [ |ll3[ |ll4| . The virtual photon can couple to any of the two nucleons, so 
that the isoscalar combinations G% m_ of th e nucleon form factors (NEMFF), 



defined in App. C, factorize to yield 115 1: 



Gc(Q') = GUQ^)xCe{u,w;Q^) 
Gq{Q') = G%{Q^)xCq{u,w;Q^) 
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Figure 8: Deuteron form factors Gc, Gq and Gm as a function of Q (updated 
from Fig. 1 of Ref. Q). The data for Gc and Gq correspond to i20 (or 
T20) mea surements (see legend of Fig. including new preliminary results 
from |109|] (triangles right). The open diamonds correspond to a second solution 
of the equations Gc,Gq = f{A, B, t2o) [fs). The Gm data correspond to the B 
measurements indicated in Fig. 0, with the addition of (full circles). The 
curves are from parameterizations I (solid line), II (dot-dashed) and III (short 
dashed) discussed in the text. 
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Figure 9: A{Q^) (solid curve) and the contributions from Gc (long dashed), 
Gq (dashed) and Gm (dotted), from parameterization I. Above 1.3 (GeV/c)^, 
the statistical significance of the data is such that the crossing of the Gc and 
Gq contributions, corresponding to ^20 = 1/V2, is not firmly established. 



Gm(Q') = GUQ^)xGL{u,w;Q^) + Gi{Q^)xCs{u,w;Q^) (36) 

The C functions are integrals of quadratic combinations of u and w. In this 
NRIA, the ratio Gq/Gc, and thus ^20, are independent of the nucleon form 
factors. The coupling of the virtual photon to the moving nucleon charges and 
to the nucleon spins both contribute the magnetization in Gm , giving rise to the 
two terms in (|3^). The elastic structure functions A{Q'^), B{Q'^) and i2o(Q^) are 



illustrated in Fig. ^ for a variety of phenomenological potentials 1 11 6 , |117 , Q , 
using the MMD parameterization of the nucleon form factors (see App. C). 
In all cases the calculations agree with one another and with the data up to 
~ 0.5 (GeV/c)^ but diverge at higher Q^. The low behaviour of the form 
factors is not determined with the same precision for each of them. Since the 
forward cross sections at low depend mostly on Gc, the slope of this form 
factor at = is determined to about 1% (see Sec. 3.1.5). The slope of Gm 



at the origin is determined by backward cross sections and is known to about 
5%. In contradistinction, the slope of Gq is known to only 15%, in the absence 
of very precise ^20 measurements at low . Still this slope is model dependent. 
To illustrate this point, we define the reduced quantity 

and show its low behaviour in Fig. At — 0, one gets 

f and ^'W(o) _ ^'9(0) rwm + ^Q^'^) rss^ 
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Figure 10: Elastic ed scattering observables for a variety of phenomenologi- 
cal potentials in the NRIA: Bonn-A (dotted curve), Bonn-B (dashed), Bonn-C 
(long-dashed), Bonn Q (dot-dashed), Reid-SC (thick dotted), Paris (solid) and 
A-vis (thick dashed). NEMFF: MMD. 
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Table 2: Dependence upon various potentials of the quantities appearing in 
Eq.(|8|). Slopes in (GeV/c)-^ 



Potential 


Gq(0) 


i20fl,(0) 


G'q(O) 


G'ciO) 


^20«(0) 


Bonn A 


24.8 


.960 


-506 


-18.1 


-1.69 


Bonn B 


25.1 


.972 


-464 


-18.1 


0.25 


Bonn C 


25.4 


.983 


-464 


-18.1 


0.47 


Bonn Q 


24.7 


.956 


-464 


-18.1 


-0.06 


Reid SC 


25.2 


.976 


-464 


-18.1 


0.32 


Paris 


25.2 


.976 


-453 


-18.5 


1.09 


A-l>i8 


24.1 


.932 


-444 


-18.1 


0.28 



Gq(0) is the model calculation of the quadrupole moment. The slope i20i?,('-') 
rather small, due to an approximate cancellation between the first two terms, 
the third one being significantly smaller. In this connection, it is instructive to 
list in Tabic |2| the model dependences of these various slopes. The value of the 
quadrupole form factor varies substantially for these calculations and is always 
smaller than the experimental value of 25.8 (p^). All of the Bonn calculations 
and the Reid-SC calculation yield the same value for the derivative of the charge 
form factor and the same of is true of the derivative of the quadrupole form factor 
with the exception of Bonn A. Therefore the spread in t'^Qj^iO) for this latter 
set comes only from the variation in the value of the quadrupole form factor. 
However, this may not be the case in general. 

The increasing disagreement between these simple calculations and the data 
at higher indicates that this approach does not contain all of the neces- 
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sary physical degrees of freedom and/or that since the momentum transfer ap- 
proaches and then surpasses the mass of the nucleon, it is necessary to consider 
the impact of special relativity on the calculation of the deuteron elastic struc- 
ture functions. We will now consider the first of these possibilities. 



5.2 Limitations of the simple potential model 



For many applications in nuclear physics, the phenomenological potential ap- 
proach is completely adequate and provides a good description of various phe- 
nomena. Its limitations were seen, however, by considering the calculation of 
electromagnetic properties of the deuteron in the preceding section. 

The inclusion of electromagnetic interactions with the two-nucleon system 
imposes the requirement that the electromagnetic current matrix elements sat- 
isfy the continuity equation 



V- < */|J(x)|^', > +^ < */|/5(x)|*, >= 



(39) 



where J(x) and p(x) are the current and charge density operators. The operators 
must then satisfy 



V-J(x)-i p(x),i7 



= 0. 



(40) 



Since either of the nuclcons can be charged, the charge density operator is 

Pi+ (41) 
while the current density takes the general form 

J = Jl + Jz + Je:r , (42) 

where pi and Jj are the charge and current densities for particle i and 3 ex is a 
possible additional contribution to the current density. 

The current of a free charged particle must satisfy the continuity equation, 
which implies 



v-J,(x)~i p,(x),r, 



Using this with (p 



(|41|) and (^) gives 

V- Je,(x)-Z [p(x),y 



= 0. 



= 0. 



(43) 
(44) 



Since pi is proportional to (1 -I- Tg ) and the two-nucleon potential V has terms 
proportional to r'^^ ■t^'^\ the second term in (^^, proportional to (r^^^ Ar*^^))3, 
does not vanish and 3ex must be nonzero. In addition, since V involves the 
coordinates of both nucleons, the current lex must also depend upon both sets 
of coordinates and is therefore a two-body operator. The physical origin of this 
contribution to the current comes from the fact that the two-nucleon potential 
contains terms corresponding to the exchange of charge between the nucleons, 
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which must in turn give rise to an associated current. These two-body currents 
are called exchange or interaction currents. 

Equation is a symmetry constraint on the theory. It cannot, however, be 
used to uniquely determine 3ex since divergenceless pieces can be added to any 
current satisfying (Q) without modifying the constraint. A unique prediction of 
the current therefore requires that the underlying dynamics of the interaction be 
specified. Viewing the nuclear force in a meson-exchange model, these two-body 
currents are naturally associated with contributions that involve the exchange of 
mesons [ |118[ . Now that it is generally accepted that quantum chromodynamics 
(QCD) is the correct description of the strong force at the scales of interest here, 
these currents must ultimately be associated with the exchange of quarks. 

This discussion leads to a consistent treatment in the context of a non- 
relativistic treatment of the two-nucleon problem. However, any attempt to 
implement this approach in a manner that is applicable to the existing data, 
for instance to the case of elastic electron-deuteron scattering, leads directly to 
consistency problems. One of the most elementary indications of this problem 
arises from the necessity of including electromagnetic form factors for the nu- 
cleons. At low momentum transfers, these form factors differ from their static 
values by terms of order Q'^/mj^. Since the dipole mass mu is of similar mag- 
nitude as the nucleon mass, this is of leading relativistic order jt? . Whenever 
the presence of the nucleon electromagnetic form factors has an appreciable 
effect on the calculations of the deuteron form factors, some effects of relativis- 
tic order are already included. It then becomes necessary to consider whether 
there are other relativistic effects of similar size that can appreciably modify the 
calculations. In fact, there are many such effects arising from Lorentz boosts 
and exchange currents that become increasingly important in the calculations 
of deuteron form factors as the four-momentum transfer increases. For this rea- 
son it is necessary to consider how relativity can be introduced consistently in 
models of the deuteron. For this purpose, we will assume that the deuteron is 
adequately described by nucleons and mesons. 

5.3 Construction of relativistic models 

The requirement of any truly relativistic model is that it must satisfy Poincare 
covariance: it must be covariant with respect to Lorentz boosts, spatial rotations 
and space-time translations. This can be imposed by requiring that the the 
operators which act as generators of these transformations satisfy the Lie algebra 
of this group: 




(45) 



{P^.P"] = 
[if\P"] = -iP^, {J\P^\ = 
\K\ P^\ = -iS''P°, [J\ PJ] = ie'^'^P' 



(46) 
(47) 
(48) 
(49) 
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where J* are the three angular momentum operators that generate rotations, 
are the generators of the three Lorentz boosts and is a four- vector contain- 
ing the energy and momentum operators that generate space-time translations. 
This differs from the corresponding algebra for the Galilean transformations in 
only two commutators. The differing commutators for the Galilean transforma- 
tions are: 

= (50) 
and 

[K'q,P^ = -iS'm (51) 

where Kq are the generators of the Galilean "boosts" and M is the mass oper- 
ator. The first of these implies that Galilean boosts commute whereas Lorentz 
boosts do not. The second is the source of the dynamical complexity of rela- 
tivistic models and theories. While the commutators for the Galilean boosts 
and the three-momentum operators are proportional to the mass, the corre- 
sponding commutators for the Poincare group are proportional to the energy 
operator, that is the Hamiltonian. Since the Hamiltonian contains the inter- 
actions between the constituents of the system, the first commutator of ( ^9|) 
implies that the Lorentz boost operators or three-momentum operators must 
also be dependent upon the interaction. 

There are two basic approaches to constructing Poincare invariant models. 
We will start with the most familiar of these, quantum field theory. 

5.3.1 Quantum field theory 

The starting point of a quantum field theory is a Lagrangian that is constructed 
to satisfy all of the required symmetries, including Poincare invariance. By 
nature field theories have an infinite number of degrees of freedom. Canonical 
quantization is performed by constructing the Hamiltonian, finding the gener- 
alized position and momentum in terms of the fields, writing the fields as an 
expansion in terms of creation and annihilation operators, and then imposing 
canonical equal time commutation (anticommutation) relations on the canonical 
variables. This yields commutation (anticommutation) relations for the creation 
and annihilation operators. An immediate consequence is that the fields com- 
mute (anticommute) for all spacelike intervals, implying that events occuring 
at spacelike separations cannot be causally connected. This is the property of 
microscopic causality or microscopic locality. It should be noted here that micro- 
causality results from imposing the commutation (anticommutation) relations 
on any spacelike hypersurface. 

The structure of field theories is very complex due to the presence of an 
infinite degrees of freedom. Since the Hamiltonian of an interacting system 
links states containing different numbers of particles, and the time-evolution 
operator of the system is given by the imaginary exponential of the interacting 
Hamiltonian, the interacting system contains contributions with any number 
of particles, from zero to infinity. Practical calculations in field theory, with 
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the exception of numerical approaches such as lattice QCD, must then intro- 
duce some method of truncating the collection of configurations contributing 
to the calculation. The most familiar approach to this problem is Feynman 
perturbation theory. Here, various configurations are carefully arranged such 
that all quantities contributing to a process, such as free propagators and ver- 
tex functions, are individually covariant with respect to noninteracting Lorentz 
transformations. The truncation then requires that there be some plausible 
scheme for systematically organizing contributions in order of their relative im- 
portance. For example, in the classic case of quantum electrodynamics where 
the coupling constant is small, the terms are ordered in powers of the coupling 
constant and truncated at some finite order. 

The drawbacks of quantum field theory for constructing models such as for 
NN interactions are related to complexities and calculational difficulties. Most 
of our knowledge of field theory is based on perturbation theory and there 
is no a priori method for determining the radius of convergence (if any) for 
a given field theory. For strong coupling theories, it is no longer plausible 
to construct a perturbation scheme in terms of simple counting of coupling 
constants and infinite sets of contributions must be summed. In practice, the 
coupling of the infinite sets of states, or n-point functions, must be truncated if 
there is to be any hope of calculation. Thus some physically reasonable scheme 
for truncation must be proposed. However, the truncation of the theory will 
usually violate some symmetries of the full theory such as crossing symmetries, 
covariance or current conservation. Local effective theories with the appropriate 
symmetries may also be nonrenormalizable. Model calculations may also include 
phenomenological elements such as form factors that are not calculated within 
the context of the field theory. This also leads to problems of consistency within 
the model and may also violate symmetries. 



5.3.2 Hamiltonian dynamics 

Although Dirac was one of the founders of quantum field theory, he soon be- 
came disillusioned with its complexity and the difhculties associated with the 
unavoidable infinities. He continued for most of the rest of his life to seek an 
alternative to quantum field theory. He assumed that the problems with field 
theory were related to starting from an unsatisfactory relativistic classical the- 
ory. He pointed to an alternate approach, starting with a theory with a fixed 
number of degrees of freedom, as is done with the nonrelativistic Schrodinger 



equation. This led to relativistic constraint dynamics |119|. In this approach, 
the dynamics of a model system is determined by choosing a mass operator 
which contains an instantaneous interaction as in the nonrelativistic potential 
theory. This basic dynamics contains a finite number of particles and has a 
corresponding Hilbert space when quantized. Covariance is imposed by con- 
structing a unitary representation of the inhomogenous Lorentz transformations 
with generators that satisfy the commutation relations of the Poincare group. 
The wave functions have the same probabilistic interpretation as in nonrelativis- 
tic quantum mechanics, but microscopic causality is not respected: the theory 
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must be constructed to respect the physical requirement of chister separabihty 
or macroscopic locahty. 

There are at least three different approaches for the quantization of such 
models. The first is the traditional method of quantizing along constant time 
surfaces (called instant form) where the evolution of the system is the usual time 
evolution which is normal to the constant-time hypersurface. The second is to 
quantize along spacelike surfaces with constant interval (called point form). The 
evolution of the system is then along a new coordinate normal to these surfaces. 
The third of these is to quantize along the light cone with evolution along the 
coordinate = x'^ + x^. Since all these choices describe surfaces with spacelike 
separations (or the infinite momentum limit of such a surface in the light-cone 
case), they are also consistent with microscopic causality for field theories, which 
may also be quantized along these surfaces. 

Particularly useful in the Hamiltonian dynamics is the fact that a careful 
construction of the mass operator can lead to equations of motion of the same 
form as the two-body Schrodinger equation. It is therefore possible to use, 
without modification, nonrelativistic potentials that have been fitted to describe 
A''A^ scattering . 

The drawbacks are related to the choice of the interaction without speci- 
fication of any underlying dynamical content. As a result, quantities such as 
electromagnetic currents can be constrained by the structure of the theory, but 
can not be uniquely determined from the interaction dynamics. 

We will now proceed to a discussion of various approaches used in construct- 
ing relativistic calculations of the elastic electromagnetic form factors for the 
deuteron. 



5.4 v/c expansions 

This approach is actually a hybrid of field theory and Hamiltonian dynamics. 
It assumes that the basic dynamical content of the deuteron is nonrelativistic 
and that the necessary relativistic effects can be described as corrections to the 



nonrelativistic current matrix elements as an expansion in v/c |12C, 121, 122]. 

The nucleon-nucleon interaction is taken to be a standard nonrelativistic po- 
tential with parameters determined by fitting to the nucleon-nucleon scattering 
data and to the deuteron binding energy. It is assumed that the potential is, 
at least in part, represented by a one-meson-exchange model, since the meson 
degrees of freedom are necessary to construct two-body exchange currents from 
simple Feynman diagrams and for constructing corrections due to retardation of 
meson propagators. Examples of the required two-body interaction currents are 
represented in the diagrams of Fig. Diagram (a) represents a contribution 
due to coupling of the photon to the current of an exchanged meson, which, 
because of G-parity, does not apply to isoscalar transitions such as ed elastic 
scattering. Diagram (b) is a so called pair current that arises due to the pro- 
jection of the interactions onto the positive-energy space only. Contributions 
that couple to the negative energy part of the nucleon propagator must then be 
included in the current. Diagram (c) is a retardation correction. Since nonrel- 
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(b) 



(c) 



(d) 



(e) 



Figure 12: Two-body currents arising in the vjc expansions (see text for expla- 
nations). 



ativistic wave functions are single-time wave functions, diagrams corresponding 
to the absorption of a photon while the meson is in flight are not included in 
the impulse approximation and must be included in the current. Retardation 
corrections can also be associated with retarded meson propagation in any of 
the other diagrams in this figure. Derivative couplings of the mesons to the 
nucleons can also give rise to contact interactions of the type shown in diagram 
(d), which is isovector and does not apply to ed elastic scattering. It is also pos- 
sible for a photon that is absorbed by an exchanged meson to excite the meson. 
An example is the p-K^ exchange current which is commonly included in calcu- 
lations of elastic electron-deuteron scattering. Such a current is represented by 
diagram (e). 

Relativity is imposed by requiring that the currents and interactions are 
consistent with operator commutators of Poincare invariance to some order in 
vjc. The dependence of the boost operators on the interaction also gives rise 
to interaction currents in addition to those characterized above. This approach 
guarantees that the interaction model can be very well constrained by data but 
its application can become technically complicated. In addition, the expansion 
in vjc must fail at some value of momentum transfer. 

The calculations shown here are from a recent paper by Arenhovel, Ritz 



and Wilbois (ARW) [123|. This paper is an extension of earlier work [116 



124 1 that included relativistic corrections for the charge and quadrupole form 
but not the magnetic. Although there is an extensive literature on 
as n icely summarized in |123|, only that paper and the earlier 



factors 
the subject 

work of Tamura [ 125 take into account all leading order terms including the 
Lorentz boost of the deuteron center of mass. The need for a realistic one- 
boson-exchange potential leads ARW to use the Bonn OBEPQ models, which 
may however not have the same precision in describing scattering data 
as more recent, but more phenomenological, potentials. Finally the Galster 
dipole parameterization for the nuclcon electromagnetic form factors was used. 
Figure ^ shows the observables ^(Q^), B{Q'^) and t2o(Q^) for this model. 
Three curves are shown for each observable, the impulse approximation (NRIA), 
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Figure 13: Elastic ed scattering observables calculated with the v/c expan- 
sion [123|: NRIA with the Bonn OBEPQ-B potential (dotted curve); same plus 
all relativistic corrections to leading order (dashed); all of the former plus the 
piTj exchange current (solid). NEMFF: Galster. 
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the impulse approximation plus all relativistic corrections to leading order, and 
all of the former plus the pir^ exchange current. The pir^ exchange current is 
purely transverse and is not required or constrained by the interaction model, 
but is the longest range isoscalar current of this type. 

The elastic structure function A{Q'^) for the impulse approximation is sub- 
stantially below the data. The inclusion of the various relativistic corrections 
brings the calculation into much better agreement with the data, while the 
addition of the pir^ exchange current has only a very small effect. 

In the case of the magnetic structure function B{Q^), the NRIA again is 
below the data and has a diffraction minimum at lower than indicated by 
the data. The addition of the relativistic corrections increases the size of the 
calculated structure function and appears to move the diffraction minimum 
above the value indicated by the data. The addition of the pir"/ exchange current 
again increases the size of the structure function to a value considerably above 
the data and presumably pushes the diffraction minimum to even higher values. 
We will return later to the question of the effect of this exchange current on the 
magnetic structure function, in the context of field-theoretical models. 

For the polarization structure function i2o(Q^), the NRIA is above the data 
for small below its minimum, and below at above the minimum. The 
addition of the relativistic corrections brings the calculation into much better 
agreement with the observed node of Gq and consequently with the data. The 
y07r7 exchange current only has a very small effect on this observable which is 
purely longitudinal. 



The examination of the figures in Ref. [123| shows that the various correc- 
tions to the NRIA are not small, they tend to be of similar magnitude, and they 
contribute to the deuteron form factors with varying signs. As a result, any 
calculation that includes some, but not all, of these corrections, is questionable. 

These calculations clearly indicate that relativistic effects, that is meson- 
exchange contributions and genuine relativistic corrections, are important even 
for relatively small momentum transfers. The domain of validity of this approach 
is thus limited, probably to < 1 (GeV/c)^, though it may lead in some 
instances to a satisfactory description of the data up to 2 (GeV/c)^ |^^. The 
study of the validity of the two-nucleon description of the deuteron down to very 
small internucleon distances, together with new data becoming available from 
Jefferson Lab, make mandatory the consideration of fully relativistic models. 

5.5 Relativistic constraint dynamics 

Figure ^ shows three recent calculations the deuteron elastic structure func- 
tions. Note that in this and following figures, the different observables are not 
shown in the same range, following the available data. This should be kept 
in mind for a meaningful comparison of models to data. The calculation of 



Allen, Klink and Polyzou (AKP) [126| uses the Argonne vis potential for an 



interaction and is quantized in point form. The calculation of Lev, Pace and 



Salme (LPS) [[7^, 127 1 uses the Nijmegen II potentia l an d is quantized in front 



form. The calculation of Forest and Schiavilla (FS) |12S] uses the Argonne vis 
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Q' (GeV/c)' 



Figure 14: Elastic ed scattering observables calculated with the Hamiltonian 
constraint dynamics. AKP |l2|, NEMFF: MMD (solid curve); LPS p^, GK85 



(dashed); FS [128|, H (dot-dashed). 



potential and is quantized in instant form. These three calculations involve 
three different approaches to construct currents. AKP and LPS use different 
procedures for constructing currents from the single-nucleon current without 
including any interaction dependent two-body currents. FS assume that the 
potential is of one-boson-exchange origin and construct the various currents as 
in the v/c expansions, but without performing this expansion. These calcu- 
lations show considerable variation in all three observables. Although the FS 
calculation works quite well for all of them, it is clear that no consensus has 
been reached concerning consistent techniques for the construction of currents 
in this framework. 
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Figure 15: Feynman diagrams respresenting the two-body scattering matrix M.. 



5.6 Field theoretical models 
5.6.1 The Bethe-Salpeter equation 

The oldest of the field-theoretical treatments of the two body-problem is due 



to Bethe and Salpeter [129|. The physical content of the Bethe-Salpeter equa- 
tion can be easily understood by considering the two-body scattering matrix 
M. ■ Figure ^ represents the Feynman expansion of the scattering matrix M 
where the nucleons are represented by solid lines and the mesons are represented 
by the dashed lines. Note that diagrams (b), (d), (e) and (f) can be reduced 
to simpler diagrams by cutting across two nucleon lines as represented by the 
dotted lines. These diagrams are said to be two-body reducible diagrams and 
the remaining ones two-body irreducible diagrams. The ability to classify all of 
the contributing diagrams as members of one or the other of these two classes 
suggests that the multiple scattering series can be resummed by separating the 
irreducible diagrams into an interaction kernel and then using an integral equa- 
tion to produce the reducible diagrams. The Feynman diagrams representing 
the two-nucleon irreducible kernel are shown in Fig. ^ The Bethe-Salpeter 
equation for the scattering matrix can then be written as 

Mip',p;P) = Vip',p,P)-i J -0^V{p',k;P)Go{k,P)Mik,p;P) (52) 

where P is the total four- momentum of the two-body system; p' , k and p are the 
final, intermediate and initial relative four-momenta of the two particles; Go is 
the free two-body propagator; and is a kernel consisting of a sum of all possible 
two-body irreducible diagrams. The Bethe-Salpeter equation is represented by 
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Figure 16: Feynman diagrams representing the two-nucleon irreducible kernel 
V. 
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Figure 17: Diagrammatic representation of the integral equation for the scat- 
tering matrix. 



the diagrams in Fig. 

Since the Feynman series is organized such that all of the individual propa- 
gators and vertices are covariant with respect to the free-particle Lorentz trans- 
formations, the sum is manifestly covariant as well. The currents can be con- 
structed by combining the free one-nucleon currents with exchange currents 
obtained by attaching a photon to every possible place within the irreducible 



Ward-Takahashi identities [130 



interaction kernel, as shown in Fig. 18. These currents will then satisfy the 



The two-body bound state is represented by the Bethe-Salpeter vertex func- 
tion O, which satisfies 

0{p, P)^^j ^(P' ^)Go(fc, P)0{k, P) . (53) 



This vertex function for the deuteron can be written | |131| as 

Ox,,baiP,P) = (r(p,P) • ^xAP)C\a ■ (54) 

P is the deuteron four-momentum, p the relative momentum of the external nu- 
cleons, (P) the polarization four- vector for the deuteron, C the Dirac charge 
conjugation matrix. The subscripts a and b are indices in the Dirac spinor space, 




Figure 18: The two-body irreducible current diagrams. 
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and can be determined by basic symmetry arguments to have the general 
form: 



r''(p,p) 




\ ^1 +171 



m 



, (55) 



where pi = ■ 
requires that 



p 

2 



+ 



p and p2 ^ Y ^ P- generahzation of the Pauli symmetry 



P'(p,P) = -CT^''^ i-p, P)C-\ 



(56) 



which places constraints on the scalar functions gi. Note that for given values 
of p and P, the vertex function depends upon eight scalar functions whereas the 
Schrodinger wave function of the deuteron is determined by the two functions 
u and w. 

The Bethe-Salpeter equation is a complete representation of all possible con- 
tributions to the two-body amplitudes for a given field theory. As such it re- 
tains all of the complexities of field theories in that there is an infinite number 
of contributions as represented by Feynman diagrams. The problem of particle 
dressing and renormalization is present, and for strong coupling theories there 
is no general scheme for organizing the renormalization program. As a practi- 
cal matter, model calculations generally assume that all propagators represent 
dressed particles with physical masses. This still leaves an infinite number of 
contributions to the irreducible kernel. Again, there is no a priori means of es- 
tablishing a reasonable truncation procedure in order to obtain a tractable set 
of contributions. However, the phenomenology of nuclear systems suggests that 
the importance of contributions to the nuclear force depends upon the range 
of the contributions due to the strong repulsion of nucleons at short distance. 
As a result, models of the deuteron using field-theoretical techniques usually 
assume that contributions can be ordered by range. In practice this reduces 
to the use of one-boson-exchange potentials. This is referred to as the ladder 
approximation to the Bethe-Salpeter equation. 

The ladder approximation violates the crossing symmetry of the two-body 
scattering amplitudes, which is a property of the full Bethe-Salpeter equation. 
This results from the elimination of the higher-order crossed-box contributions 
to the kernel. A related problem is that the two-body equation no longer re- 
duces to a one-body wave equation at the limit of infinite mass for one of the 
constituents (defined as the static limit) [^32j. 

A further difhculty for producing models of the deuteron using the field- 
theoretical methods is that phenomenology often requires the introduction of 
nonrenormalizable couplings such as the derivative coupling of the pion and 
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tensor coupling of vector mesons. The loop integration of the Bethe-Salpeter 
equation must then be cut off to eliminate unphysically large or infinite contri- 
butions. This is usually done by introducing form factors for strong coupling 
vertices, which also attempts to include finite-size effects for the hadrons. In 
addition, the hadrons have a finite electromagnetic size so that electromagnetic 
form factors are also required by the phenomenology. The introduction of form 
factors can, however, result in violations of gauge invariance. 

The solution of the Bethe-Salpeter equation is also a calculational challenge 
since it is a four-dimensional integral equation with a complicated analytical 
structure. It has however been solved for the two-nucleon system in Euclidean 
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space 133 and applied to the calculation of electron-deuteron scattering [134 



5.6.2 Quasipotential Equations 

The Bethe-Salpeter equation is a four-dimensional integral equation. As such 
it is much more difficult to solve numerically than the comparable nonrelativis- 
tic three-dimensional Lippmann-Schwinger equation. To simplify the solution 
of the relativ istic equations , one reso r ts t o the infinite class of quasipotential 
equations ||, |l36[ [S^, |l3|, |l3|, |l4^, fil], [l4| |l4|. These equations are re- 
lated to the Bethe-Salpeter equation (|52|) by replacing the free propagator Gq 
by a new propagator g. The scattering matrix equation can now be formally 
rewritten as 

M{p',p;P)^U{p\p;P)-i J -0^U{p' ,k; P)g{k; P)M{k,p; P) , (57) 
where U is the quasipotential defined as 

Uip',p;P) = Vip',p;P)-i J -^Vip' ,k-P){G^{k-P) ~ g{k-P))U{k,p-P) . 

(58) 

This pair of equations is equivalent to (52) and represents a re-summation of the 



multiple scattering series. The new propagator g is usually chosen to include 
a constraint in the form of a delta function involving either the relative energy 
or time in such a way as to reduce (^7|) to three dimensions. In addition g 
must be chosen such that it has the same residue as Go along the righthand 
elastic cut. This guarantees that the discontinuity of U produces only inelastic 
contributions as in the case of (^2|). 

Although it appears that the reduction of ( pTj) to three dimensions is of great 
practical advantage, it should be noted that (|58| ) is still a four-dimensional 
integral equation of difficulty comparable to (|5^. The real utility of these 
equations comes about when ( |58|) is approximated by iteration and truncation. 
Then U can be obtained by quadrature of ( p8| ) and used in the three-dimensional 
equation (|5^). 

It would appear that we have achieved a reduction in the difficulty of solution 
of the problem at the expense of introducing considerable additional ambiguity 
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since there is an infinite number of possible choices for g which satisfy the above 
specified requirements. However, this freedom is turned to an advantage when 
noting that g can be chosen to maximize the convergence of (|58|). In essence, 
the parts of the ladder and crossed box diagrams which tend to cancel are being 
summed in the quasipotential equation ( |58| ) rather than in the equation for the 
scattering matrix ([57|). Indeed, the lowest order calculations using a variety of 



quasipotential prescriptions have been calculated for scalar particles |144| and 
shown to be closer to the results of the complete Bethe-Salpeter equation than 
is the result from the ladder approximation to the Bethe-Salpeter equation. In 
addition, it has been shown that there exists an infinite number of quasipotential 
equations which have the correct static limit. This is also a refiection of the 



improved convergence of these quasipotential equations | 145 | 



The results from three different quasipotential models are presented here. 



The first of these is the calculation of Hummel and Tjon (HT) 1 146 ] based on the 



BSLT equation 136 . In this case, the propagator gsSLT contains a term d{po) 
which forces the relative energy of the two nucleons to zero. One can say that 
the two nucleons are equally off mass-shell. The vertex functions are calculated 
using the BSLT equation with a one-boson exchange interaction. They still 
have eight components. The current matrix elements are then calculated by 
replacing the Bethe-Salpeter vertex functions with the BSLT vertex function, 
assuming that it is energy independent. In the examples shown here the nucleon 
electromagnetic form factors are taken to be the Hohler 8.2 parameterization. 



The second model is that of Phillips, Wallace and Devine (PWD) [147 



148, 14£] where a one-boson-exchange interaction is used with the single-time 
equation that introduces a constraint that the relative time be zero. This is very 
close to the HT approach since this choice means that the propagator gET is 
independent of the relative energy of the two nucleons. The single-time deuteron 
vertex functions still have eight components. A consistent treatment of the 
current is constructed to guarantee current conservation, but Lorentz covariance 
is violated. The MMD parameterization of the nucleon electromagnetic form 
factor are used. 



The third model is that of Van Ordcn, Devine and Gross (VODG) |150| , |51 



152 1, based on the Gross equation [141, 142 [. One of the nucleons is placed on 



its positive energy mass shell, which gives: 



9Gross OC 

Since this constraint is itself manifestly covariant, the Gross equation ampli- 
tudes are manifestly covariant. But it is not symmetric in the nucleons and the 
exchange symmetry must be recovered by symmetrizing the interaction kernel. 
This procedure introduces unphysical singularities, which are treated in princi- 
pal value to avoid unitarity problems and have little numerical effect in a weakly 
bound system such as the deuteron. As a result of placing one nucleon on shell 
in (|5^), the Gross equation deuteron vertex function has four-components that 
can be represented in terms of an S wave, a D wave and singlet and triplet P 
wave functions. 
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A systematic procedure for constructing the effective current operators for 
the Gross equation exists |153[ and has been shown to be rigorous to all orders 



and remarkably robust under trucation |154]. The result is that Ward identities 



for the Gross equation are maintained and the calculation is gauge invariant. 
A one-boson-exchange interaction is used, that has been fit to NN scatter- 



ing data up to a lab kinetic energy of about 300 MeV |155, |156[ . The fit is 
reasonable, but not at the level of precision obtained with modern fully phe- 
nomenological potentials. 

Several unique features appear in this model. All of the strong vertices are 
multiplied by a product of three form factors, 

h{p'')h{p^)f{{p' -p)'). (60) 

Here p and p' are the initial and final four-momenta of the nucleon and p' — p 
is the four-momentum of the meson. The meson form factor is taken to be 

/(^') - .. . , (61) 



and the nucleon form factor 



' (A2 - p2)2 + (A2 _ j„2)2 ■ 1° ) 



The prescription of Gross and Riska |153| is used to construct a single- 
nucleon electromagnetic current which is consistent with the strong- vertex form 
factors and phenomenological singic-nuclcon electromagnetic form factors. This 
current is constructed in such a manner that the one- and two-body Ward 
identities are of the same form as for the local field theory. The simplest form 
that this single-nucleon current can take is 

J«^(/,p) = F,{Q^)Up'\p^)^^ + ^^h,{p'\p^)^a'^'^q, 

^FMQW.P')^-,'-^. (63) 

where the nucleon form factors are 

F.{Q') = liFfiQ') + F^iQ^)-T,) . (64) 
The factors that depend upon the nucleon momenta are 

JO[P ,p )- ^^^,2^ p2_p/2 + /j(p2) p/2_p2' ^O^^ 

and 

,2 2.fh{p') h{p'^)\ 4m2 
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Figure 19: Elastic ed scattering observables calculated with Quasipotential 
equations. HT-RIA, NEMFF: H (dotted curve); PDW-RIA, NEMFF: MMD 
(dashed); VODG-CIA, NEMFF: MMD (solid). 
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The function /io(p'^,p^) must be equal to one when both nucleons are on mass 
shell, but is otherwise unconstrained. For simplicity, the calculations shown here 
use ho — fo- The form factor FioiQ"^) must obey -Fio(O) = 1, but is otherwise 
unconstrained and is chosen to be of a dipole form. 

Figure ^ shows the elastic structure functions for these models calculated 
in the relativistic impulse approximation (RIA) . As expected from the nature of 
the quasipotential approximations, the calculations of HT and PWD lead to very 
similar results for all three observables (in spite of the fact that they use different 
NEMFF). The VODG calculation uses the MMD form factors. For A{Q'^), it 
is systematically larger than the others due to additional contributions that 
are necessary to conserve the current within the context of the Gross equation. 
These additions yield to the so-called Complete Impulse Approximation (CIA). 

For B{Q^), the HT and PWD calculations are systematically below the 
data and have a minimum at a lower than is indicated by the data. The 
VODG calculation is higher and has a minimum at higher than the other 
two calculations. This difference has been shown to be due to a small P wave 
component of relativistic origin that contributes to the normalization of the 



vertex function at the level of 0.01% |152 . This clearly indicates the sensitivity 



of the position of the Gm node to small components. 

Finally, in the case of t2o(Q^); the three models have a very similar behavior. 

The sensitivity of the observables to the choice of single-nucleon form factors 
is shown in Fig. |2^, using the VODG calculation in the CIA approximation. All 
but the V02 parameterization are standard form factors that appear in the 
literature. The V02 form factor is modeled after the MMD form factor but 
adjusted to fit the new data on G^/G^^ from Jefferson Lab [157|. As expected 
from the discussion in App. C, A{Q'^) and B{Q^) are very sensitive to the choice 
of single nucleon form factor while t2o (Q^) is almost totally insensitive to it. The 
various j4(Q^) curves in Fig. ^ may be directly related to the corresponding 
G% curves in Fig. Note that for the relativistic approaches, in contrast to 
the NRIA, the nucleon form factors cannot simply be factorized, so that there is 
no reason that i2o(Q^) be completely independent of the NEMFF. Much more 
accurate data on the single-nucleon form factors that will be forthcoming from 
experimental facilities around the world will, hopefully, provide much greater 
constraints on these form factors. 

Figure ^ shows the effect of the addition of p-n^ exchange currents to the 
three calculations. HT were the first to implement the calculation of the pir^ 
contribution to the deuteron form factors in a relativistic model. They also 
included a oja^ exchange current, which compensates their large effect of the 
pTT^. This last contribution is not present in the calculation presented here, 
because there is considerable uncertainty about the a meson, and consequently 
about its couplings and vertex form factors (see 1 158 for a recent discussion of 
the (7cij(T7 coupling constant). Two sets of calculations for the VODG structure 
functions arc shown, one with the MMD form factors as before, and one with 
the V02 form factors. Comparing to Fig. |l^, it is clear that this exchange 
current can result in considerable variation in A{Q'^). The variation in B{Q^) 
is much smaller, in part because of the smaller range considered. The pnj 
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Figure 20: Deuteron elastic structure functions calculated using the CIA of 
VODG with various parameterizations of the nucleon form factors: Galster (dot- 
ted curve), Hohler (long dashed), GK85 (dot-dashed), GK92 (dashed), MMD 
(solid), V02 (thick dotted). 
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exchange current shifts the node of Gm to lower values of for all calculations, 
in contrast with the nonrelativistic calculations using the v/c expansion where 
this current has the opposite effect. This is the result of the tensor coupling 
of the p to the nucleon, which appears to be a higher order contribution in 
the v/c expansion, but weighted by a large coefficient. In this way, expansion 
schemes may lead to wrong estimates without a careful examination of higher- 
order contributions. Finally, the addition of the pir^ contribution has a small 
effect on i2o(Q^), but brings the calculations in a somewhat better agreement 
with the available data. 



The variation in the pir^f contributions from the various models |l4q, 159 



160 , |161| is associated with ambiguities in these contributions. Although the 
coupling constant for the pi:^ vertex is constrained by the decay width T{p ^ 
TTj), little is known about the fall off of the associated form factor. Figure ^ 
shows a number of different form factors for this vertex. HT use the VMD 
form factor which is the hardest of the available form factors while VODG use 
the Rome 2 form factor which is the softest one. PWD use an intermediate 
parameterization. The A{Q'^) data favor the use of the softest possible form 
factor. In contradistinction, a recent theoretical reexamination of the pirj form 



factor [162 1 results in a dependence close to the VMD one. Whether the pirj 
contribution to the deuteron form factors is really small, or suppressed by other 
contributions is still an open question. 

Another source of ambiguity, explicitly present in the the VODG calculation, 
is the choice of the form factor Fio(Q^) in (p3|). In these calculations, i^io(Q^) 
is adjusted to optimize the fit of the calculation to the data for the deuteron 
structure functions. By using a very hard form factor it is possible to obtain 
an extremely good fit to the data. This along with ambiguities in the pirj 
form factor means that no absolute predictions of the structure functions can be 
obtained unless some means can be found to physically constrain the ambiguities 
in the models. 

5.6.3 Light-front field theory 

The above approaches based on the Bethe-Salpeter equation are generally dis- 
cussed in the context of Feynman perturbation theory with equal time quanti- 
zation. Microcausality, however, only requires that the theory be quantized on 
a spacelike hypersurface. A special limiting case of such a surface is the light 
cone where the spacetime interval approaches zero. Quantizing field theories 
on the light cone has long been common practice in describing deep inelastic 
scattering where the large four-momenta kinematically favor contributions to 
scattering very near the light cone. Typically, the light-cone approach is or- 
ganized as a "time-ordered" perturbation expansion in terms of the light-front 
time x" = x'^ ~ x^. The wave functions are then Fock-space wave functions 
with a probabilistic interpretation as in the case of Schrodinger wave functions. 

A particular problem with light-front field theory is that the conventional 
choice of a fixed light-front orientation violates manifest covariance, although 



matrix elements remain covariant. One solution to this problem |163 is to de 



45 



10" 



10" 



10"' 



10" 



10" 



Tfe 1 ' 


1 1 1 1 1 

A 


1 ' ; 


10"^ 
















10"= 








10"*= 
















10"^ 








10"= 


1 


1,1, 


1 





1 2 

0.5 
0.0 
-0.5 
-1.0 
-1.5 



3 4 5 6 




- 1 1 1 1 1 1 


1 1 1 1 1 1 1 




- \ji 

, Cf"Pr 1 ^ , 


f 

1 , , , 


1 , , , - 


0.5 


1.0 

Q' (GeV/c)' 


1.5 



Figure 21: Elastic ed scattering observables calculated with Quasipotential 
equations, including the pn^ contribution. HT (dotted curve), PDW (dashed), 
VODG-MMD (solid), V0DG-S02 with i^m off-shell form factor adjusted (thick 
dotted). 
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Figure 22: Dif ferent prescriptions for t he pn'y form factor: VMD (dott ed c urve), 
Gross-Ito llll (dashed), Rome 1 (long dashed), Mitchell-Tandy |6l| (dot- 
dashed), Rome 2 |l6§ (solid). 



scribe the quantization surface in terms of an arbitrary direction on the light 
cone given by a light-like four-vector uj'^ and require that the quantization sur- 
face be described by a; • a; = 0. This renders the theory manifestly covariant, at 
the expense of the additional complication that the wave functions and opera- 
tors become dependent on the vector us. This is covariant light-front dynamics 
(LFD). Matrix elements and observables, in a complete calculation, do not de- 
pend on Lo. In a practical calculation, a definite procedure to eliminate the non 
physical w-dependent terms is applied. 



The deuteron elastic structure functions from this approach |163, 164, 165| 
are shown in Fig. 

In these calculations, the deuteron wave function has six components (fi), 
two of which (/i and /2) correspond to the usual S and D components in the 
nonrelativistic limit. These components are calculated perturbatively from the 



usual nonrelativistic wave function. Schematically |166| 



^(k,n;/i,/2,/3,/4,/5,/6)~ Jv(k',k,n,Md)^PNRO^;u,w)d''k' , (67) 

where k is the relative momentum of the two nucleons and n a unit vector along 
the three-vector uj. The NN interaction kernel is of a one-boson-exchange 
type, calculated within the framework of LFD using the parameters (such as 
coupling constants) of the Bonn potential. The single-nucleon electromagnetic 
form factors are the MMD parameterization (the effect of various NEMFF was 
considered in [ |165| ). This calculation produces a reasonable description of the 
data for A{Q'^) and i2o(Q^) but the minimum in B{Q^) occurs well below the 
position indicated by the data. Given the sensitivity of the magnetic form factor 
to small effects, this may simply be the result of the perturbative treatment of 
the small components of the wave functions in this calculation. 
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Indeed, work is in progress to improve the calculation of the interaction 



kernel V and of the complete wave function |167 . Meson exchange contributions 
are only partially included in the original calculations. The addition of recoil 
terms and of the pnj contribution was also initiated: the recoil terms move 
the node of Gm to higher Q^, while the pn^ contribution does not modify 
appreciably any observable up to ~ 2 (GeV/c)^ [168|. A complete calculation 
using LFD seems at hand and very promising. 

5.6.4 Effective field theory 

In principle, all these field theoretical techniques should yield the same results, 
given the same dynamical input. This statement applies to the calculated matrix 
elements which are the physical observables of the field theory. The identifica- 
tion of wave functions and operators varies from formulation to formulation 
and do not have unique meanings. Furthermore, within any given formulation, 
unitary transformations of the Lagrangian, such as field redefinitions, can also 
move contributions between the wave functions and operators. Such ambigu- 
ities are already present in nonrelativistic models, which may be shown to be 



equivalent, at least in part, up to a unitary transformation [ [44[ 16E |. None of 
this would be of particular concern if it were not necessary to truncate all of the 
approaches for reasons of practicality. Consequently, the physical content of the 
various formulations varies, resulting in a variety of ambiguities, some of which 
have been discussed above. A major problem then is that there is, in general, 
no organizing principle present in these calculations which indicates the relative 
importance of various physical contributions to guide the choice of truncation 
schemes. 

A promising development that could help to resolve this problem is the 



application of effective field theories to nuclear systems [17C, 171 



172, 173, 174 



175, 176, 177, 178|. The basic idea of effective field theory is that at low energies 
the observables of a theory are largely insensitive to the details of short-range 
contributions to the interactions. The long-range degrees of freedom are then 
treated in detail and the short range pieces are replaced by contact interactions. 
The Lagrangian is written as a sum of terms containing contact interactions with 
increasing numbers of derivatives of the fields. This constitutes an expansion of 
the theory in terms of momenta which are small compared to some scale chosen 
to separate the short and long range physics. Contributions to observables 
are then ordered according to this small momentum, providing a well defined 
counting scheme that controls approximations to the theory. 

For nuclear systems, the appropriate effective field theory is Chiral Pertur- 
bation Theory (xPT) since the long-range part of the nuclear force is associated 
with the pion, a Goldstone boson. A considerable amount of effort is being ex- 
pended in the application of xPT to low energy nucleon-nucleon scattering and 
to the deuteron. The NN system is particularly challenging in that its scat- 
tering lengths are large, although its bound state, the deuteron, is only weakly 
bound. This implies that there is a dynamical scale in the problem in addition 
to the chiral scale. The application of xPT to the deuteron has followed two 
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approaches. The first [170| appHes the chiral counting scheme at the level of 



the interaction kernel. This approach appears to converge, but is cutoff depen- 



dent at each order. The second |173, 174| applies the counting scheme at the 



level of the scattering matrix and current matrix elements and treats the pions 
perturbatively while iterating the contact interactions. It gives results that are 



cutoff independent, but does not appear to be converging [179|. The problem 
of maintaining a consistent counting scheme in the infinite sums of diagrams 
necessary to describe the bound states has not yet been resolved. From the 
standpoint of the various models used at higher momentum transfers, it is to be 
hoped that effective field theory will provide some insight into the organization 
of the various approaches. 

5.7 Deuteron models with nucleon isobar contributions 

In all of the approaches and models discussed to this point, the assumption 
has been made that the only relevant degrees of freedom are the nucleons and 
mesons. However, excitations of the nucleons to isobar states may produce con- 
tributions to the interactions of the same range as the heavier mesons. The 
deuteron wave function is in this case modified to include, in addition to the S 
and D wave NN components, AA and A'^A'^*componcnts. A neutrino experi- 
ment, though subject to interpretation, indicates an upper limit of 0.4 % to the 



amount of AA components in the deuteron |18C]. 

Examples of two calculations including these additional isobaric components 
are shown here. In both cases, the basic model involves interactions due to one 
meson exchange which can couple the nucleon-nucleon channel with channels 
containing isobars. The quark model is used to relate isobar-meson couplings 
(for instance AAtt) to the corresponding nucleon-meson couplings. The fully 
coupled (nonrelativistic) system is then solved for the deuteron bound state 
and the nucleon-nucleon scattering states. Figure shows the calculations of 



Dymarz and Khanna (DK) |181] and Blunden et al. (BGL) [182|. In both cases 
presented here, only the A(1232) isobar is included, and the A electromagnetic 
form factors are assumed to be proportional to the NEMFF. The DK calculation 
contains AA components '^Si,^ DiJ DiJ Gi with a total probability of 0.36%. 
The elastic form factors include contributions from single-nucleon currents with 
the GK form factors, nucleon pair and pnj exchange currents and the isobar 
current contributions. Two models are shown for BGL, both with only ^Di 
and ^ Di A A components.. Two different quark models are assumed to fix at a 
given radius ro the boundary conditions used in the determination of the wave 
functions. Model C' has ro = 0.74 fm and yields an isobar contribution of 1.8%, 
while model D' uses the Cloudy Bag Model to fix ro at 1.05 fm and results in 
an isobar contribution of 7.2%. Meson exchange contributions are included as 
well. 

Within these models, t2o is very sensitive to the amount of A A components 
and favors the smallest probability. Although the BGL model C' with the Hohlcr 
NEMFF seems to give an adequate description of the data, there is considerable 
ambiguity in these calculations due to a lack of knowledge of the isobar magnetic 
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moments and form factors along with the difficulty of completely constraining 
the isobar contributions from the NN scattering data. It is clear, however, that 
isobar components must ultimately be considered in any complete calculation 
of deuteron electromagnetic properties in the context of hadronic models. 

5.8 Quarks and gluons 

5.8.1 Nonrelativistic quark models 

Ultimately, many of the ambiguities associated with the hadronic model calcu- 
lations can only be removed by describing the deuteron in terms of the funda- 
mental quark and gluon degrees of freedom. In the absence of the capability to 
directly solve the QCD Lagrangian at low energies for the deuteron, we are left 
to explore possible QCD effects in the context of quark models. This greatly 
increases the difficulty of treating the NN system since what is a two-body 
problem in the hadronic models, becomes at least a six-body problem for quark 
models. Calculations for two examples of quark cluster models are presented 
here. These calculations are very similar in concept and are based on a quark 
cluster model using the resonating group method to describe the interaction 
of the two three-quark clusters. The quarks interact via a quadratic confining 
potential and a one-gluon-exchange potential. Long range interactions are also 
provided via tt and p exchange between quarks. The calculation of Buchmann, 



Yamauchi and Faessler (BYF) [183| contains only the tt exchange interaction. 



while that of Ito and Kisslingcr (IK) |184| contains both meson-exchange con 



tributions. This approach naturally contains currents associated with the in- 
dividual clusters as well as exchange currents associated with quark exchange 
between clusters. Since the cluster wave functions are derived using oscillator- 
like confining forces, the cluster form factors tend to have a Gaussian form. This 
is clearly not consistent with the data for singlc-nucleon form factors and in both 
cases the electromagnetic cluster form factors are replaced by phenomenological 
NEMFF. The results are shown in Fig. |2^. Neither of these calculations provides 
an adequate description of the data, but given their necessary simplicity, they 
are remarkably close to the data. In the BYF calculation, the impulse approx- 
imation is not as reliable as in NN models because of the simple modeling of 
the intermediate range interaction, but the size of the genuine quark exchange 
contributions, due to the antisymmetrization of the six-quarks wave function, 
is significant enough to affect a comparison with the data. 

5.8.2 Perturbative QCD 

This subject was reviewed recently jl2). The success and shortcoming of pertur- 
bative quantum-chromodynamics (PQCD) applied to the asymptotic behaviour 
of form factors may be better illustrated now with recent data from Jefferson 
Lab. We will only recall briefly the main predictions: 



Dimensional scaling |185|: this property can be derived within QCD, but 



is more general and was established before this theory. It is based on 
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Figure 24 
components 



Elastic ed scattering observables for models including nucleon isobar 
DK, 



NEMFF: GK [|l81[ (solid curve); BGL-C, NEMFF: H |82| 
(dotted); BGL-C, NEMFF: GK (dashed); BGL-D', NEMFF: GK (dot-dashed). 



52 




q' (GeV/c)' 



Figure 25: Elastic ed observables calculated with quark-cluster models. IK [184| 
(dotted curve); BYF [|l8|: lA (dashed); BYF: lA + MFC (dot-dashed); BYF: 
former -t- quark exchange contributions (solid). 
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the hypothesis that the momentum transfer is shared among ah the con- 
stituent quarks of the system. In the case of a hadron composed of n 
quarks, one would expect the leading form factor to behave asymptoti- 
cally as Q"^^""^). Qualitatively, this can be viewed as the probability 
of having n — 1 quarks within a transverse distance of 1/Q of the quark 
struck by the virtual photon. Since this must be true both in the initial 
and final states, one expects ~ Q"'*^""^) and for the deuteron (n = 6), 

Logarithmic corrections to the leading amplitude, together with the fac- 



torization of the nucleon form factors in the weak binding limit, yield |186|: 



where Fn is the nucleon form factor, A an energy scale characteristic of 
QCD (A ~ 200 MeV) and e a calculable number much smaller than 1. An 
attempt to directly compute A and its normalization in PQCD was not 



successfull [12, 187 



• Helicity conservation at each photon/gluon-quark vertex implies that the 
dominant contribution to elastic electron deuteron scattering should come 
from the configuration where the deuteron has helicity zero in both ini- 



tial and final states [188 . Expressing the form factors in an helicity basis 
in the light-cone frame, this is equivalent to the prediction that G^Jq, the 
component of the current matrix ele ment between states of helicity 



0, is the leading amplitude. It was argued [189| that this helicity conserv- 
ing amplitude should dominate the ed scattering for ^ 2Amd — 0.8 
(GeV/c)'^. Other components G^q (single helicity flip) and Gj_ (double 
helicity flip) should be suppressed by respectively one and two powers of 
Q: 



G+o = a ( - ) G^ and G|_ = b[-] G^^ when Q ^ oo (69) 



A\ „^ . „^ , /A 

Q 

Since the usual form factors are linear combinations of the G^^, 's, the 
asymptotic behaviour of observable ratios such as B/A, t2o and t2i can 
easily be calculated. Note that the logarithmic corrections such as appear- 
ing in ( |6^ ) may only be calculated for the dominant, helicity conserving, 
form factor. Equation (^9|) thus assumes the same logarithmic corrections 
for all helicity amplitudes. 

How do these predictions compare with recent ed data ? The A data, now 
extending up to = 6 (GeV/c)^ [^, is suggestive of the expected Q"^" 
behaviour, though a Q^^^ behaviour is still not excluded (see Fig. ^): a fit 
using the five highest data points in Rcf. [^ to the dependence A ^ Q^^™ 
yields m = 8.0 ± 0.6. Excluding the less precise last point yields m — 8.7 ± 0.7. 
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Figure 26: A{Q'^) in a Log-Log scale. The data is from JLab/HallA (9§. The 
straight hnes illustrate a Q"^® (dashed) and a Q^^" (long-dashed) dependence. 
The solid line represents Eq.(p8|). All curves are normalized to the point at 
Q2 = 4.95 (GeV/c)2. 



More interestingly, using the dipole form factor for Fjy in (|68|), the behaviour 
of A is reproduced between 2 and 6 (GeV/c)^. As for the helicity amplitudes, 
their behaviour is tested only up to ~ 2 (GeV/c)^, which is the range of the 
available B 



96 1 and t2 



dat a. The pure dominance of the helicity 0^0 
transition (a = 6 = in (|69|) 189 ) does not account for this data (see Fig. p7|) . 
The prescription a — 5 and 6 = |190 was built to generate a node in Gm, but 
in order to qualitatively reproduce the t2j data as well, one is led to a = 1.8 
and 6 = 38 |191|. This implies that the double helicity flip amplitude is as large 
as the non helicity flip amplitude and contradicts the applicability of PQCD in 
the momentum range considered. Fits to data using these helicity amplitudes 

All these simple ansatz lead to a sharp 
2.5 (GeV/c)^, which may be an interesting 



192 



lead to the same conclusion | 178| , 
increase of the ratio B/A for 
feature for planned experiments at JLab. 



> 



5.9 Further comparison between models and data 

To conclude this section, we recapitulate some of the theoretical results, in com- 
parison with observables (Fig. |2^ ) and data on separated form factors (Fig. p9|) . 
In the case of A and i?, deviations with respect to an average representation of 
the data (parameterization I) are presented. A summary of remaining ambigu- 
ities and foreseable progress will be given in Sec. ||. 
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Figure 27: Ratios B/A, <20 and <2i from light-cone frame helicity amplitudes. 
The curves correspond to Eq .(|69|) with different values of a and b (see text), 
from Refs. |l89| (soUd Hnc), [l90| (dashed) and ||l9l| (dot-dashed). Note the 
logarithmic scale in Q^. 
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Figure 28: AA/A and AB/B, in %, and ^20 as a function of Q^. The deviations 
of A and i? are calculated with respect to parameterization I. A data legend: 
Saclay @ (open diamonds), SLAC Q (open circles), JLab/HallA Q (filled 
circles), JLab/HallC (filled diamonds). B data legend: see Fig. ^. ^20 
data legend: see FigsT^ andj. Theoretical calculations: NRIA + MFC + 
relativistic corrections (Sec. ^) |123| (d ashed curve) and (d otted), CIA 
+ MEC (Sec. ^I^ ) updated froirTlLSlt (sohd), LFD (Sec.^sTe^ |l65| (dot- 
dashed). 
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6 The nucleon momentum distribution in the 
deuteron 



Several processes have been investigated for which the cross sections, in the 
nonrelativistic impulse approximation, are proportional to the distribution of 
the internal momentum p in the deuteron, 

p[p)^u^{p)+w''{p), (70) 

with u{p) and w(p) defined by Eq.(^). Note that p, the conjugate variable 
of the relative coordinate r, is half the relative momentum between the two 
nucleons, so that an experiment probing a momentum p may be related to elastic 
electron scattering at a momentum transfer q = 2p. These processes include 
quasi-free scattering on one of the nucleons, d{p,2p)n, d{e,e'p)n and d(e, e'), 
or equivalently the detection of the spectator proton in high energy deuteron 
hadronic break-up, Aid^p) . Indeed, they all are proportional to p{p) up to 
p ~ 200 MeV/c , but deviations from the impulse approximation occur for higher 
values of p . Final state interaction, dynamical excitation of the A resonance 
and in inclusive reactions pion production have then to be taken account, thus 
rendering the interpretation of the experiments less straightforward. 

Another interesting feature of these processes is the possibility to access the 
ratio w{p)/u{p) with the measurement of deuteron tensor polarization observ- 
ables. A simple relationship between the analyzing power T20 and the ratio 
w{p)/u{p) can be derived in the (nonrelativistic) impulse approximation, and 
again deviations are seen, or are to be expected, above p ~ 200 MeV/c . 



6.1 d(e,e'p)n measurements at high missing momenta 

The study of the single-particle properties of nuclei through (e, e'p) reactions 
are the subject of an excellent review 0|. Concerning the deuteron, the (neu- 
tron) missing momentum is identified, in the plane-wave impulse approximation 
(PWIA), with the internal momentum p . The experiments reaching the highest 
missing momenta have been carried out at Saclay (500 MeV/c) ]193| |, NIKHEF 
(700 MeV/c) pi and MAMI (950 MeV/c) |jl9|l, but in kinematical conditions 



which are not always optimal to study the high momentum components in the 
deuteron wave function. Deviations from the PWIA are of the order of 50 % at 
500 MeV/c and can reach a factor 10 at 1 GeV/c (see Fig. ^0|). They can be 
understood, if only qualitatively at the highest missing momenta, in terms of 
final state interactions, meson exchange currents and A excita tion. 

Other d{e,e'p)n measureme nts (see the reviews [ 196 , 197 and, among the 
most recent experiments, Ref. |198|) do not address specifically the subject of 
the nucleon momentum distribution in the deuteron. Their theoretical under- 
standing is however crucial for a comprehension of this subject. 

The first (e, e'p) measurements using a tensor polarized target have been 
carried out at Novosibirsk and NIKHEF [pOO| . They were limited to 

low values of missing momenta because of the available luminosity and of the 
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Figure 30: Ratio of the experimental d{e, e'p)n to the calculated P WIA (using 
the Paris potential) cross sections, corresponding to the data of Ref. [195|. Sys- 
tematic errors are of the order of 10 %. Note the change of vertical scale at 
p = 400 MeV/c . 



detector acceptance. This situation will improve with the planned experiments 
using the BLAST detector in the Bates stretcher ring [ p01| . In the PWIA, these 
experiments are directly sensitive to the ratio w{jj)/u{jj) (see below Eg. ([71 



6.2 d(e,e') measurements and y-scaling 

The concept of nuclear y-scaling gives a significant insight in the determination 
of momentum distribution in nuclei |202| . In the inclusive electron scattering off 
deuterium, for sufhciently large values of momentum transfer and in the PWIA, 
the cross sections are expected to be proportional to the electron-nucleon cross 
sections and the proportionality factor /(y) = 27r Jj^i pp{p)dp represents the 
longitudinal (along q) momentum distribution. Properly taking into account 
final state interactions (a correction of about a factor 2 above 300 McV/c), 
a momentum distribution was extracted in a nonrelativistic formalism up to 



p = 600 MeV/c 203 . It is remarkably close to the one calculated from the Paris 



potential, but the authors caution about the absence of relativistic corrections. 



6.3 Hadronic deuteron break-up at high energy 

Several reactions using high energy (polarized) deuteron beams should allow a 
study of the dpn vertex, provided the deuteron dissociates "cleanly", that is 
with only one nucleon participating in the process. In the exclusive reactions 
p{d,p)d, p{d,pp)n as well as in the inclusive A[d^p)X, in some kinematical 
conditions, the dominant process is the one where the fast forward proton (or 
one of the detected protons) can be treated as a spectator. The energies of the 
deuteron beams available at SATURNE and at Dubna (several GeV), together 
with the relatively high hadronic cross sections, have provided a handle on these 
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processes involving high momentum components in the deuteron wave function. 
Moreover, the tensor analyzing power of these reactions, in the nonrelativistic 
impulse approximation, depends only on the D/S ratio at a given momentum: 

w{w~2^/2u) 
T20 oc — — . (71) 



Note the similarity between Eqs.(^^ and (71). The former is a function of 
277Gq/3Gc, while the latter depends on w/u, resulting (in the impulse ap- 
proximation) in similar shapes for T20 in ed elastic scattering and in deuteron 
break-up. 

A relativistic treatment using the deuteron in an infinite momentum frame 



allows to define a new variable k [204 used as the argument of u and w instead 
of p for the calculation of the observables. In other words, the argument of the 
wave function is no longer equal to the momentum of the spectator nucleon. 
Once this transformation is performed, all inclusive data A{d,p)X (see the re- 



view [205 1 and Refs. [20(:, 207|), for different beam energies and target nuclei, 
scale approximatively as a momentum distribution p{k). T20 however deviates 
significantly from expectations within the impulse approximation as of fc ~ 200 
MeV/c. This fact has sometimes been interpretated as the signature of non 
conventional components in the deuteron wave function, but this interpretation 
is not compatible with the behaviour of i20 in elastic ed scattering. In reality, 
final state interactions and pion production alter significantly the interpreta- 



tion of these experiments (see e.g. [208|). The polarization transfer from vector 



polarized deuterons to the fast protons ha s als o been measur ed up to k ~ 600 
MeV/c. The exclusive channels p{d, _p)d §, |o|| and p{d, pp)n [|09[ |lC | are like- 



wise difficult to interpretate unambiguously for internal momenta larger than 
300 MeV/c. 

7 The deuteron as a source of "free" neutrons 

The neutron being loosely bound in the deuteron, deuterium has often been 
used as a substitute for a neutron target and deuteron beams as a source of 
neutron beams. We briefly mention the list of these applications : 

1. ed elastic scattering to extract the neutron charge form factor 0"% , 

2. Quasi-elastic ~ct {it , e' n)p or dllf , e'lt)p to measure G^/G'lj , 

3. Quasi-elastic d{e,e'n)p and d{e,e'p)n to measure G\.jlG\j , 

4. Deep inelastic scattering of leptons to extract the neutron structure func- 
tions. 

5. At intermediate to high energies, the break-up of a vector polarized deuteron 
beam can be used to obtain a beam of polarized neutrons. A high inten- 
sity may be achieved using the inclusive reaction ~S -(-^ Be — > + A". For 
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a better definition of the neutron energy, the beam may be tagged by the 
detection of the spectator proton: p{a ,pTt)p. 



In all cases but the first one 211 1, corrections due to the deuteron structure 



and to the reaction mechanism were shown to be either negligible or reliably 
calculable (see [1135| for example as an application of relativistic techniques 
discussed in Sec. |5.6.l|) . When a polarized target is used (cases 2 and 4) or when 
using polarized deuteron beams (case 5), the effective neutron polarization is 
equal to the deuteron vector polarization multiplied by (1 — 3Pd/2) . This factor 
accounts for the fact that, because of the deuteron D state, the neutron spin is 
not always aligned the deuteron spin. 



8 Prospects for the future 

The study of the electromagnetic properties and form factors of the deuteron, 
from the birth of nuclear physics to the advent of hadronic physics, taking 
into account the internal structure of the nuclcons and their excitations, has 
been very rich. Yet it is not completed. Ambiguities have been and still are 
pointed out along this path: although the descriptions of processes involved 
can be satisfactory, the calculations and observables do not in general allow 
for a unique determination of the NN interaction, of the neutron charge form 
factor, of the isoscalar meson exchange currents and of the precise dynamics 
of the system. As already pointed out, some characteristics like the off-shell 
behaviour of the NN interaction are not strictly speaking observables, and as 
such cannot be determined unambiguously. As for the manifestation of the 
underlying quark substructure in the nuclear properties, it is as elusive as ever. 
However the progress in experiment and theory, as summarized in this paper, 
is impressive and allows for promising perspectives: 

• The nucleon electromagnetic form factors are being measured with a re- 
newed precision. In particular, the poorly known neutron electric form 
factor is being determined at various laboratories with polarization 
techniques which make this measurement independent of the deuteron 
structure. In a few years, all four NEMFF will be better known up to 

~ 2 (GeV/c)^. The description of the the deuteron form factors to 
higher four-momentum transfers is now affected by measurements of G^. 
Once these are completed, a small remaining ambiguity due to the poor 
knowledge of G^ above 2 (GeV/c)^ will still be present. Yet new param- 
eterizations of the NEMFF, guided by theoretical models, should become 
available soon and be taken into account in future calculations of the 
deuteron form factors. 

• Modern NN interaction models are now fitted directly to the NN elastic 
scattering data and reach a high degree of precision. But, when compared 
to some models of the 1980's, these are more empirical. The lack of knowl- 
edge (or assumption) of the underlying dynamics prevents the calculation 
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of such effects as the Lorentz boost of the dcuteron wave function. It 
would be highly desirable to have potentials of the one-boson exchange 
type brought to the degree of precision of the modern phenomenological 
potentials. This is also true for OBE potentials used in completely rcl- 
ativistic calculations. In addition, these potentials should reproduce the 
elastic NN scattering data up to 600 MeV of kinetic laboratory energy in 
order to match the Q^-range of the existing ed elastic scattering data. 

There is a definite procedure to construct rolativistic corrections, starting 
from a nonrelativistic model. Still, most calculations in the past have 
neglected one correction or another. 

From a theoretical point of view, the most satisfying success in the past 
twelve years is the implementation of various fully relativistic calcula- 
tions. Quasi-potentials approximations to the Bethe-Salpeter equation 
have been applied with success to the calculation of the electromagnetic 
deuteron form factors. The comparative validity of each of these approx- 
imations should be studied further. Likewise, light-front dynamics yields 
results which compare very favourably with data and will still be im- 
proved. Concerning Hamiltonian constraint dynamics, calculations using 
the three different forms of quantization invoked by Dirac have been de- 
veloped. The success of QPE/BSE and LED has reaching consequences 
beyond the structure of the deuteron. The relevance and applicability of 
relativity in nuclei can now be explored in the A=3 systems. The same 
techniques are also applied to qq configurations. 

The pirj isoscalar meson exchange contribution, and possibly the ua^ or 
other shorter range processes, still remain difficult to evaluate reliably, 
because of the lack of constraints on the associated form factors. Within 
the nucleon-meson picture of the deuteron, ed elastic scattering may pro- 
vide a way to determine these form factors, provided all points above are 
addressed in a systematic way. 

The role of nuclcon isobaric excitations is still very much model dependent. 
Still, within the existing models, the ed elastic scattering data does not 
favour very sizeable AA components in the deuteron wave function. 

The recent ed elastic scattering data from Jefferson Lab, reaching now the 
highest possible four-momentum transfers for the ^20 and A observables 
(respectively 1.7 and 6 (GeV/c)^), are still compatible with the descrip- 
tion of the deuteron in terms of nucleons and mesons. This is somewhat 
surprising since internucleonic scales of 0.1-0.4 fm are being probed. These 
distances are smaller than the size of the nucleons themselves, and pre- 
sumably of the same order of magnitude as the nucleon quark cores. Still 
no distinctive experimental or theoretical signature of the manifestation 
of quarks in this process was identified. 
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• Quarks degrees of freedom may be explicitly taken into account at inter- 
mediate four-momentum transfers via models, and at high via pertur- 
bative QCD. If the recent j4(Q^) data seem to follow the Q^-dependence 
expected from PQCD, it is not so for B{Q^) and t2o(Q^), albeit at lower 
Q^. An absolute determination of the leading PQCD amplitude would 
clearly be desirable but this depends on a reliable description of the soft 
parts of the amplitude which are not calculable in PQCD. As for quark 
models of the deuteron, they seem to indicate a specific role played by 
quark exchange processes between the nucleons at short distances. Un- 
fortunately, these models are nonrelativistic by nature and the predicted 
effects occur at a scale where relativity should be taken into account. Due 
to the explicit appearance of the quark-gluon degrees of freedom, this a 
much more difficult problem than for the meson-nucleon models. Appre- 
ciable progress in quark models of the deuteron will require substantial 
improvements in the technology of relativistic many-body physics. 

• The planned measurements of t2o at Bates and Novosibirsk will locate 
more precisely the position of the node of the charge monopole form factor, 
while the ones of B{Q^) at Jefferson Lab will be performed around and 
beyond the node of the magnetic dipole form factor. 

• The recent focus on intermediate and high should not be detrimental 
to the required precision in accounting for the low data and static 
properties. An updated experimental status of the latter was given. The 
low behaviour of any given calculation should be checked carefully, for 
example in such representations as in Figs |l^ and 

Going beyond elastic ed scattering, the same systematic expansions in v/c, or 
fully relativistic models, will be applied to the calculation of the electromagnetic 
form factors of the A=3 nuclei, and of the deuteron electrodisintegration. This 
will provide additional contraints on the remaining ambiguities inherent in the 
meson-nucleon models and will lead to a more coherent understanding of the 
relativistic structure of few-body nuclei. 
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A Beyond one photon exchange at high ? 

Since the ed elastic cross sections are now measured up to very high momentum 
transfers, the question of the validity of the one-photon exchange approximation 
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Figure 31: Two-photon exchange diagram. 



(used throughout this review) has to be examined. As shown in Sec. 4.4, the 
set of ed elastic scattering data, as parameterized up to ~ 2 (GeV/c)^, is 
compatible with Eqs. (25-31), thus providing evidence for the vaUdity of the 
one-photon exchange approximation, at least within experimental errors. A 
two-photon exchange contribution would affect the definition of all observables 
presented in Sec. 4.2, introducing additional kinematical factors and structure 
functions [|| 



Several papers have addressed this question in the early 1970's [212, 213 



214 |215|| : the two-photon exchange process may be viewed in Glauber theory 
as the double scattering illustrated in Fig. Estimates of this amplitude 
were made with the assumption that each of the exchanged photons carries 
half of the momentum transfer. With such approximations, and depending 
on the relative phase of the two amplitudes, the double scattering amplitudes 
could change the ed forward cross sections by as much as 10% for ~ 1 — 2 
(GeV/c)^. A complete calculation would involve an integration over all possible 
four-momenta of the exchanged photons and over all possible intermediate pn 
states. This has not been attempted yet, except for small and for a spinless 
deuteron [216|. In this particular case, definite corrections to the deuteron 
charge radius could be made. Finally, the high experimental data on the 
structure function A ^ Q were recently compared, but no conclusive 
evidence for a two-photon exchange contribution could be reached |217|. New 
dedicated experiments and calculations would certainly clarify this important 
question. 



B Polarized deuteron targets - Polarimeters 

The parallel development of polarized atomic beam sources and solid polarized 
targets originated in the late 1950's. The use of adiabatic transitions between 
hyperfine energy levels of atomic hydrogen or deuterium for the former and the 
dynamic polarization of nuclei for the latter were decisive improvements due in 



good part to the work of A. Abragam [218|. Their use in electron scattering 
experiments is however much more recent. 

The polarized sources led to the acceleration of polarized deuteron beams. 
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Figure 32: The polarimeter POLDER |221| as used in the Jefferson Lab t2o 
experiment |10S]. 



which in turn were used to conceive and cahbrate deuteron polarimeters. There 
are however very few exemples of efhcient deuteron tensor polarimeters at inter- 
mediate energies (in ed elastic scattering, the recoil deuteron energy is related 
to the four-momentum transfer by = Q"^ /2Md)- The most recent ones were 



developed at the synchrotron SATURNE-2 [gl9|: the AHEAD polarimeter [g20| 



was based on dp elastic and inelastic scattering in the 100-200 MeV range, while 



POLDER (Fig. |32| and |221]), using the charge exchange reaction p((i,pp)ri, was 
operated between 160 and 520 MeV. There is at present no good concept to 
develop a tensor polarimeter for deuteron energies above 500 MeV. In con- 



tradistinction, vector polarimeters are available [219 , but they are less useful 
for the determination of the deuteron form factors (see Sec. |4.2| and A 
peculiar feature of the double scattering experiments using a tensor polarime- 
ter is related to the nature of the ^20 moment : contrarily to other moments 
and to the more familiar vector polarization, it does not induce any azimuthal 
dependence in the distribution of counts in the polarimeter, which goes like 

Noi0P°') X [1 + t2QT^°\0P°^) + 2t2iT|°'(6'P°') cos(pP°' + 2t22Tl°\eP°^) cos2(pP°'] . 

One then needs to know the absolute response of the polarimeter to both po- 
larized and unpolarized deuterons in order to extract t20- This requires the 
separate polarimeter calibration with a deuteron beam of known polarization. 
In some instances, the angular dependence of T'I'q if large enough, may be used 
to avoid or check the determination of an absolute normalization ^ |103[ . 

Atomic beams have also been used as internal targets in storage rings. In 
this case the thinness of the target is compensated by the multiple passage of 
the circulating electron beam. The target thickness was later increased by a 



factor of about 100 by accumulating the polarized atomic beam in a cell |222| 



Alternative plans to use an internal target based on the spin-exchange optical 
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pumping technique [223| were developed at Argonne and Novosibirsk, but not 



implemented 224 



Cryogenic solid polarized targets (ND3 or deuterated buthanol) cannot with- 
stand much electron beam intensity, but this has significantly improved in 



the past few years 225 . The target tensor polarization or alignment {P^z 
or A) is generally rather small: it is related to the vector polarization by 
Pzz = 2 — ^4 — 3P^ I lie I and the obtained values of Pz are at best of the 



order of 0.5. 

All these techniques were applied to ed elastic scattering for the measure- 
ment of t2o or T20. They are compared quantitatively in Table ^, using a figure 
of merit defined for polarimeter experiments as F = E'^ LVLe.leiT2Q^Y ^-nd for 
polarized target experiments as = E'^LVt^A^ /2. All quantities are defined in 
Table ^. The factor E"^ accounts for the approximate energy dependence of the 
ed cross section for a fixed value of Q. We have also indicated in that table our 
estimation of the presently achievable (but not planned) highest figure of merit 
with polarized targets, either internal in the upgraded HERMES ]226| | configura- 
tion at HERA, or external at Jefferson Laboratory. With present day technology 
and ideas, it seems very difficult to extend the tensor polarization measurements 
beyond the four- momentum transfer values reached in the JLab/POLDER ex- 



periment |103|. But upcoming experiments |109, 201| will cover an intermediate 



Q-range with different systematic uncertainties than in the past polarimeter 



experiments ||77|, 103], and in some cases better statistical precision. 



C Nucleon electromagnetic form factors 

We present here the most commonly used parameterizations of the nucleon form 
factors (NEMFF) in various deuteron form factors calculations. In this case the 
quantities of interest are the isoscalar combinations G% = j^+G^ j^. We 
recall that in the NRIA, the deuteron charge form factors are proportional to 



G% ; in this case, the A^ elastic structure function (Sec. 4.4) is proportional to 



(Cf;)^ while t2Q is independent of any nucleon form factor. 
The dipole parameterization 

The proton form factors and the neutron magnetic form factor can be written 
approximatively as 

GI = ^ = ^=Gd={i + 78^) ' (72) 



Mp Mn V 0-71 

while the neutron electric form factor is taken to be zero or following the Galster 
parameterization |88[ : 

Gl - -T^^^nGD (73) 

where a and b are two free parameters adjusted on data. Frequently used values 
of these parameters are given in Table 0. 
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Recoil polarimeters 

ed —^ e~S 


Polarized targets e~S ed 


Internal 


External 




Bates/ 


Bates/ 


JLab/ 


VEPP-3 


NIKHEF 


Bates/ 


Achievable 


Bonn 


Achievable 




Argonne 


AHEAD 


POLDER 






BLAST 










[76] 


[77] 


[103] 


[105]/[109] 


[108] 


[201] 




[110] 




Beam energy E (GeV) 


.371 


.850 


4.05 


2.0 


.704 


1.0 


27.5 


2.0 


6.0 


Beam intensity 


30 ^lk 


30 AiA 


100 AiA 


200/100 mA 


150 mA 


100 mA 


30 mA 


0.4 nA 


80 nA 


Target thickness (d^m'^) 


4 X 10^1 


2.4 X 10^^ 


6 X 10^^ 


6 X 10"/10^'^ 


2 X 10^^ 


10'^ 


10'^ 


9 X 10^^ 


2x10^'-^ 


Luminosity L (cm"*s~^) 


8 X lO^'^ 


4.5 X 10^^ 


4 X 10^** 


7 X lO^'VS X 10^1 


2 X 10^1 


5x10^^ 


2x10^^ 


2 X 10^^ 




Solid angle He (msr) 


20 


18 


6 


120/160 


150 


100 


150 


5 


6 


Loss factor I 


.5 


.25 


.5 














Polarimeter efficiency e 


10-4 


2 X 10-^ 


2 X 10-^ 














Analyzing power TI'q ' 


-0.8 


-0.4 


-0.2 














Target polarization 








.46/.85 


.83 


.85 


.85 


.12 


.21 


Figure of merit F 


3 X 10^« 


5 X 10^1 


1.4 X 10^^ 


8 X 10^**/2 X 10^1 


9 X 10^** 


dxlff'^ 


1.5xl(fi'-^ 


6 X 10^« 




Qmax (fm^i) 


2.0 


4.6 


6.7 


2.9/4.0 


3.2 


(4.5) 


(7) 


(3.6) 


(6.5) 



□ 



□ 



□ 



Table 4: Neutron electric form factor parameters for Eg. ([73|) . 



a 


b 


Rcf. 


1. 


5.6 




8S 




1.25 


18.3 
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1. 


3.4 


[22' 


1 



An alternative parameterization of Gd is a product of two different monopole 
form factors |228 
form factors [|65|. 



and G^j can also be described as a sum of four monopole 



The IJL-G parameterization 

The IJL parameterization [ p29[ rests on the vector meson dominance (VMD) 
model: the interaction of the photon with the nucleon is mediated by the p 
(isovector) and by the ui and (p (isoscalar) mesons. The parameters are adjusted 
on data available in 1972 and the form factors do not have the asymptotic be- 
haviour anticipated from PQCD, so that this parameterization should not be 
used beyond 1 (GeV/c)^. In several calculations of the deuteron form factors, 
a parameterization using IJL for G^ iG^/ G^j, and the Galster parameter- 
ization ( [t^ ) for G^ was often used. It is denoted by IJL-G. 



The Hohler parameterization 

This parameterization |230| also rests on VMD, but uses a larger and corrected 
data set, as well as additional heavier mesons. 



The GK parameterization 



In additi on to the VMD model (with the p and uj mesons) , this parameterization 
(GK85) |231] incorporates scaling la ws at high momentum transfer compatible 
with PQCD. A later version (GK92) |232|] includ es t he 4> meson but is in worse 
agreement with recent precise data on G^^/G^^ [157|. 



The MMD parameterization 



This more recent parameterization [233 includes four isovector and three isoscalar 
mesons and offers in addition a more precise description of the low Q behaviour 
of the form factors. 



Very recent data (see for example the reviews in [234[) on the proton and 
neutron form factors have not yet been incorporated into any of these parame- 
terizations, which differ appreciably from each other. At the present time, the 
MMD parameterization seems the most satisfactory, but more precise data and 
parameterizations of the NEMFF are certainly needed 
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are 



Finally, in Fig. 

plotted the isoscalar nucleon form factors, which are the relevant combinations 
for the deuteron form factors calculations. 
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Figure 33: Ratios of the isoscalar nucleoli form factors (electric and magnetic) 
to the dipole form factor. The parameterizations are the dipole-Galster of 
Eqs.(|7|j7|) (dotted hne), IJL-G (dashed), Hohler (long dashed), GK85 (dot- 
dashed) and MMD (solid). See text for notations and references. 
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